Pairing and Density Correlations of Stripe Electrons 
in a Two-Dimensional Antiferromagnet 
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We study a one-dimensional electron liquid embedded in a 2D antiferromagnetic insulator, and 
coupled to it via a weak antiferromagnetic spin exchange interaction. We argue that this model may 
qualitatively capture the physics of a single charge stripe in the cuprates on length- and time scales 
shorter than those set by its fluctuation dynamics. Using a local mean-field approach we identify 
the low-energy effective theory that describes the electronic spin sector of the stripe as that of a 
sine-Gordon model. We determine its phases via a perturbative renormalization group analysis. For 
realistic values of the model parameters we obtain a phase characterized by enhanced spin density 
and composite charge density wave correlations, coexisting with subleading triplet and composite 
singlet pairing correlations. This result is shown to be independent of the spatial orientation of the 
stripe on the square lattice. Slow transverse fluctuations of the stripes tend to suppress the density 
correlations, thus promoting the pairing instabilities. The largest amplitudes for the composite 
instabilities appear when the stripe forms an antiphase domain wall in the antiferromagnet. For 
twisted spin alignments the amplitudes decrease and leave room for a new type of composite pairing 
correlation, breaking parity but preserving time reversal symmetry. 

PACS numbers: 71.27.+a, 71.10.Hf, 74.20.Mn 
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I. INTRODUCTION 



Extensive experimental studies - including elastic and 
inelastic neutron scattering^ ARPESS, /xSR£, and NMR 
experiments^ - have confirmed that stripe formation is 
a property common to most high-T c cuprates. In the 
underdoped regime, at some critical hole doping, the 
mobile holes segregate into an array of "stripes" that 
slice the copper-oxide planes into alternating phase - an- 
tiphase antiferromagnetic domains. The stripes coexist 
with superconductivity, but as one enters the overdoped 
region they begin to evaporate, signaling a crossover to a 
conventional metal with a uniform charge distribution. 
Significantly, stripe phases are observed also in other 
doped antiferromagnets, such as the "Nickelates"— and 
the "Manganites"— (colossal magnetoresistance materi- 
als, where the "stripes" are actually 2D sheets of hole- 
rich regions). This suggests that stripe formation is a 
robust and generic property of this class of matter. Still, 
the basic questions why stripes form and what role they 
play for superconductivity in the cuprates remain con- 
troversial. 

Early mean-field calculations on the 2D Hubbard 
modefi suggested that the stripe phase is due to the re- 
duction in kinetic energy of holes propagating transverse 
to the stripes. In this approach, however, the possi- 
ble connection to superconductivity is left unanswered. 
In an alternative approach 8 , it is argued that stripes 



form as a response to the competition between long-range 
Coulomb interactions (which push the holes apart) and 
short-range antiferromagnetic interactions (which tend to 
"phase separate" the holes into a single region). Within 
this scenario it has been argued that a proposed spin gap 
from the undoped domains is transmitted to the stripes 
via pair hopping ("spin gap proximity effect"—), lead- 
ing to enhanced charge density wave (CDW) as well as 
superconducting pairing correlations along the stripes. 
For static stripes (as seen in e.g. in the Nickelates, or 
the Nd-doped L^-zSr^NdCuGi) the CDW correlations 
dominate. In the presence of transverse stripe fluctua- 
tions, however, these appear to die ouliiS, possibly open- 
ing a door to superconductivity. Other scenarios, where 
the stripes actually compete with superconductivity, have 
also been proposed^. 

Most of the theoretical attempts to explore the prop- 
erties of stripes model these as a collection of one- 
dimensional (ID) or quasi-lD electron liquids 12 , cou- 
pled to their neighbors 13 , or to an insulating back- 
ground, either via pair hopping of charge carriers (as 
in the "spin gap proximity effect"—) or by a spin ex- 
change. The various spin-exchange scenarios that have 
been suggestecU^i&iiSi also predict that a spin gap opens 
in the spectrum of the stripe electrons, signaling en- 
hanced superconducting fluctuations along the stripes. 
In fact, it is common to find a dynamically generated spin 
gap for a one-dimensional electron gas (1DEG) coupled 
to an active environmenliSiiLiSiiiL 2 ^, of which an antifer- 
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romagnet is a particular realizationisii 5 ^^!* 2 ^ 2 ^^ 2 *. 

The simplest such model is maybe that of the ID 
Kondo-Heisenberg lattice (KHL) which consists of a 
1DEG interacting weakly with an antiferromagnetic 
Heisenberg spin-1/2 chain by a Kondo coupling. Away 
from half-filling this model has a spin gap^SiS and one 
thus expects the presence of superconducting correla- 
tions. Indeed, it was shown recently 25 - 2 — that the spin 
gap supports composita^^ odd-frequency odd-parity 
singlet pairing 29 as well as a composite CDW. A general- 
ization of this model that may mimic stripe physics more 
closely is that of a 1DEG coupled by a Kondo coupling to 
two non-interacting antiferromagnetic Heisenberg spin- 
1/2 chains, together emulating the insulating background 
in which a stripe is embedded. Rather surprisingly, as 
shown recently 3 ^, this generalized model has no spin gap 
but instead renormalizes to a fixed point belonging to 
the class of chirally stabilized electron liquids^!. Still, 
the model exhibits the same unconventional pairing in- 
stabilities as found for the ID KHL 25 i 26 . 

A "strong" interpretation of the results in Refs. 25,26 
and 30 may seem to exclude spin exchange as a possi- 
ble source of the spin gap in the high-T c cuprates: odd- 
frequency pairing appears difficult to reconcile with the 
experimental observation that superconductivity in these 
compounds is due to d-wave BCS paired electrons. How- 
ever, the recent report that underdoped BSCCO breaks 
time reversal symmetry - in the "normal" as well as the 
superconducting state - cautions us that the case may 
not be closed. The time reversal breaking is seen below a 
temperature T gap at which a pseudogapS opens, suggest- 
ing that it is connected with some order parameter that 
develops enhanced correlations below this characteristic 
temperature^. It has been argued that the pseudogap in 
the cuprates may be identified with the amplitude of the 
pairing order parameter, with long-range superconduct- 
ing order appearing at the onset of global phase coher- 
ence (carried by Josephson tunneling of pairs between 
the stripes)— One may envision a variant of this sce- 
nario where spin exchange between the stripes and their 
environment (maybe in conjunction with pair hopping) 
supports two coexisting types of quasi-one dimensional 
pairing correlations below T gap , one of which breaks time 
reversal. As one approaches the superconducting transi- 
tion, the enhanced stripe fluctuations may favor the other 
type (which could re-emerge as long-range d-wave order 
via the dimensional crossover— implied at T c ), while the 
channel that exhibits time-reversal breaking remains in- 
coherent, with only finite-range correlations present. Al- 
though speculative only, the viability of this brand of 
scenario can be judged only by more closely examining 
the physics driven by a stripe-environment spin-exchange 
interaction. This is the purpose of our paper. 

We shall consider an extended version of the model in 
Ref. 16, where a ID electron liquid (representing a sin- 
gle stripe) is embedded in a 2D antiferromagnetic back- 
ground, and coupled to it via an antiferromagnetic spin 
exchange. We show that this setup leads to a spin-gap 



phase for the electrons on the stripe, and we identify its 
leading instabilities. We further address the question to 
what extent the instabilities found are sensitive to the 
relative orientation of the staggered magnetizations on 
each side of the stripe. In the simplest case of a site- 
centered stripe^ the spin alignment across the stripe is 
antiferromagnetic (phase- antiphase domains). However, 
the alignment is not expected to be perfect and it is 
therefore important to check the stability of the spin gap 
phase w.r.t. deviations from the phase- antiphase orienta- 
tion of the magnetic domains separated by the stripe. In 
addition, we shall explore the issue whether the spatial 
orientation of the stripe on the underlying lattice may 
influence the stripe electron dynamics when the domi- 
nant interaction with the environment is that of a spin 
exchange. 

This latter question is of particular relevance con- 
sidering recent experimental findings of "diagonal 
stripes" in the underdoped glassy phase of the cuprates. 
As discovered by Wakimoto et aZ— , the insulating 
Lai.95Sro.o5Cu04 exhibits sharp two-dimensional elastic 
magnetic peaks at (ir ± e, n ± e) (in tetragonal square 
lattice notation, with e ~ x ~ 0.05, x being the doping 
level). Assuming that the magnetic peaks are associated 
with charge stripe order—, this implies that static stripes 
run along the diagonal of the square Cu 2+ lattices that 
make up the CuO-planes in this compound. This is in 
exact analogy to the diagonal static stripe structure seen 
(and theoretically predicted^) in the insulating Nickelate 
La2- x Sr :E Ni04+ 2 ;, but different from the structure in 
superconducting I^-xSro.osCuO^ (with x > 0.05) 
where the stripes are oriented along the copper-oxide 
bonds ("collinear stripes"). Very recently, these finding 
were extended to the full insulating spin glass phase in 
La 2 -xSr x Cu0 4 (0.02 < x < 0.055)22. Thus, the onset of 
superconductivity in the low-temperature underdoped 
region appears to proceed via a stripe rotation by 45°, 
from a diagonal to a collinear stripe configuration. In 
the case of a collinear (site-centered) structure, a stripe 
is embedded in a local antiferromagnetic background. In 
contrast, in the diagonal structure the stripe electrons 
experience a local ferromagnetic environment. In both 
cases the antiferromagnetic domains are 7r-shifted across 
the stripes, and it is a priori not clear whether the 
different local orderings along the stripes influence the 
electron dynamics differently To clarify the situation 
requires a careful study, and we here make a first 
attempt on it. 

To isolate the core of the problem we shall make a few 
simplifying assumptions: 

(i) We study the electron dynamics on a single stripe, 
and, in the first part of our analysis, neglect its possible 
coupling to neighboring stripes. Moreover, the stripe is 
taken to be static. This implies that for a fluctuating 
stripe (as typically seen in a superconducting phase) we 
can only hope to capture processes on length- and time 
scales shorter than those set by its fluctuation dynamics. 
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This is expected to be much slower than the dynamics of 
charge carriers along the stripe: The latter appear on an 
energy scale ~ leV, whereas the stripe fluctuations are 
coordinated with those of the localized spins, at a scale 
~ 1- lOmeV. Having obtained the characteristic features 
of a single static stripe we then add "by hand" the trans- 
verse fluctuations and interstripe couplings, and study 
their effect on the pairing- and density correlations of 
the stripe electrons. 

(ii) We model the stripe - at Fermi momenta incommen- 
surate with the underlying lattice - as a one- dimensional 
metallic wire. Thus, we assume that the disorder (from 
e.g. dopant potentials) is sufficiently weak so that local- 
ization effects set in on length scales much larger than 
those that we probe here. 

(iii) There is an important distinction between site- 
centered and bond-centered stripes^l: The spin align- 
ment across an anti-phase domain wall is antiferromag- 
netic for a site-centered stripe and ferromagnetic for a 
bond-centered stripe (which has a finite width). Moti- 
vated by recent theoretical work^fi where bond-centered 
configurations appear to be inconsistent with ARPES 
studies of La2_ a: Sr a ;Cu0^i, we here focus on the simplest 
case of a site-centered charge stripe, to be described by a 
ID Hubbard model. We shall explore elsewhere the case 
of bond-centered stripes, building on the corresponding 
analysis by Krotov, Lee and Balatsky«i£> of a Hubbard 
ladder in an antiferromagnetic environment. 

(iv) As suggested by neutron scattering data on the rel- 
evant materials^, the environment is Neel ordered up to 
some characteristic scale (which, in the relevant temper- 
ature range, is much larger than the linear dimension of 
a stripe), with a 7r-shift across the stripe when this is 
site centered (phase - antiphase domains). In our for- 
mal analysis we depict each Neel ordered domain as a 
semi- infinite 2D Heisenberg antiferromagnet, and ignore 
possible topological effects that may be present for finite- 
width insulating domains, or "spin ladders"—. We shall 
give precise estimates for the range of validity of this ap- 
proximation, thus establishing its physical relevance. 

(v) As we have already discussed, we couple the 
stripe electrons to its insulating environment exclusively 
through a spin exchange interaction. Given that the 
Fermi momentum of the stripe is incommensurate with 
that of any low-lying excitation of the environment, ex- 
cursions of single charge carriers is a process that violates 
momentum conservation, and hence is suppressed (on the 
time scales defined in (i)). Pair hopping is still allowed, 
provided that the pair carries zero total momentum. As 
suggested by the analysis in Ref. 9, pair hopping is fa- 
vored as a dominant process when the low-lying spin ex- 
citations of the environment are gapped^. When such 
a gap is absent, as is the case when the environment is 
Neel ordered, the virtual hybridization between delocal- 
ized levels on the stripe and the localized levels in the 
environment produces an effective spin exchange that is 
expected to compete effectively with pair hopping. Here, 
we focus on the effect of the spin exchange. 



(vi) We limit our attention to the physically relevant 
case of a weak spin exchange Jk between stripe and en- 
vironment, < Jk -C Jh, where Jh is the antiferromag- 
netic exchange between the localized spins in the environ- 
ment. This allows us to treat the problem in a continuum 
limitfi. 

(vii) Finally, we stress that finite-size or boundary ef- 
fects of the ID electrons^ are not included in our anal- 
ysis. As the stripes in the cuprates are mesoscopic 
structures^, these effects should in principle be taken 
into account. However, as they are not expected to qual- 
itatively change the conclusions arrived at in the large- 
distance limit considered here, we leave this study for the 
future. 

Clearly, by assumptions (i) - (vii) we loose several 
facets of the full problem. Still, we believe that our 
"stripped-down" approach has its merits: Not only does 
it isolate and expose a crucial element of " stripe physics" , 
but as we shall show, it allows us to perform a well- 
controlled analytical study, producing results that can be 
taken as a reliable starting point for more realistic stud- 
ies. Moreover, the problem as defined by (i)-(vii) is im- 
portant in its own right, and is of relevance to the more 
general issue of one-dimensional electron liquids in active 
environments^i2*iS*i2iSii. This is a central problem in the 
theory of correlated electrons, motivated by experiments 
on quasi-lD organic conductors^, nanowires^i, and edge 
states in quantum Hall systems^. 

Given the assumptions (i) - (vii), we model the stripe 
by an extended (U — V) Hubbard chain, weakly cou- 
pled to a phase - antiphase antiferromagnetic environ- 
ment by a Kondo lattice interaction. Treating the Hub- 
bard chain via standard bosonization while describing the 
environment by a nonlinear er-model, we follow the ap- 
proach introduced in Ref. 16 and exploit the symmetry 
breaking in the magnetic environment to "absorb" the 
Kondo lattice interaction as an effective spin-spin inter- 
action among the stripe electrons. In this way we obtain 
an effective low-energy model for the stripe electrons - 
decoupled from the environment - and accessible to a 
well-controlled perturbative renormalization group anal- 
ysis. This allows us to pinpoint the dynamic instabilities 
in the low-energy, weak-coupling (Jr- <C Jh < \U\ , |V|) 
limit. 

Our most important results can be summarized: 

• For realistic values of the model parameters, and 
with a phase-antiphase Neel-configuration across 
the stripe, an electronic spin gap opens on the 
stripe with a spin-density and a composite charge- 
density wave as the leading instabilities. The sub- 
leading instability is that of conventional triplet 
pairing, coexisting with composite singlet pairing 
(which breaks parity and time reversal). Slow 
transverse stripe fluctuations tend to suppress the 
density correlations, thus promoting the pairing in- 
stabilities. 

• The low-energy physics is insensitive to the spatial 



orientation of the stripe on the lattice: The results 
summarized above hold for both collinear and di- 
agonal stripes (with the possible exception that the 
composite singlet pairing is suppressed for a diag- 
onal stripe). 

• The instabilities found for the phase-antiphase 
Neel-configuration are still present when the rel- 
ative orientation of the staggered magnetizations 
on the respective sides of a (collinear) stripe has 
been twisted by an arbitrary angle. In addition, 
the twist allows for a novel type of composite pair- 
ing correlations to appear, respecting time reversal 
but breaking parity. 

The paper is organized as follows: In Sec. II we intro- 
duce the lattice models for site-centered collinear and di- 
agonal stripes coupled to a phase-antiphase environment 
by a Kondo interaction, and derive the corresponding 
low-energy effective actions. In Sec. Ill we perform an 
RG analysis and identify the order parameter correlations 
along the stripes that get enhanced by the spin exchange. 
In Sec IV we then study - for the case of a collinear stripe 
- the stability of these correlations w.r.t. perturbations 
of the relative orientation of the spin alignments across 
the stripe. Sec. IV, finally, contains a summary and a 
brief discussion of our results. Throughout the paper we 
try to supply sufficient information to make the analysis 
accessible also to the non-expert. 

II. THE MODEL 

For clarity, we shall treat the collinear and diagonal 
stripe configurations separately. The effective low-energy 
theories that emerge in the two cases are essentially the 
same (with certain provisos), but to arrive at this re- 
sult requires some care. Much of the analysis builds 
upon well-known results, but in order to make the exposi- 
tion self-contained we outline the most important points. 
Also, some key elements are new or need particular at- 
tention. 



A. Collinear Stripes 

We represent the stripe (running, say, along the x- 
direction of a square lattice) by an extended (U — V) 
Hubbard chain, H Hubbard, coupled via a Kondo lattice 
interaction H-Kondo to the nearest localized spins on each 
side of the stripe. These spins, like the rest of the lo- 
calized spins, interact mutually via an antiferromagnetic 
nearest- neighbor Heisenberg spin exchange Hafm, and 
reside in one of the two semi-infinite antiferromagnetic 
domains that surround the stripe, denoted A and B re- 
spectively (see FIG. 1). 

The A and B domains are assumed to be antiferromag- 
netically ordered, and correlated via a 7r-shift across the 
stripe (phase - antiphase domains) but there is no direct 
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FIG. 1: Collinear stripe structure. 



interaction between A and B spins^. Thus, we study 
the lattice model 

H = H Hubbard + ^ ( H AFM + H Kondo)' i 1 ) 
i=A,B 

where 

Huubbard = ~t ^ ( c l+l,a C r,a + h.C.) 
r, a 

+ U n r ^n r ^ + V h r n r+ i , (2) 



H-afm — Jh 2_^i (S r jd) ■ S r+1 j(i) + S r j(i) ■ S r j(i) +1 ), 

J H >0, i = A,B, (3) 

Hxondo — Jk ^ c\, a a a f3C r ^-{S r ,j(A) =l +S rt j(B) =l ), 

Jk>0. (4) 

Here c r . Q is a stripe electron operator at site r with spin 
index a =t,|, f>>r,a = c r a c r,a is the number operator, 
n r = J2 a ™r,a is the density operator, and Syjra is the 
operator for a localized spin at a lattice site with coordi- 
nate r (jW > 1) in a direction parallel with (transverse 
to) the stripe. The vector of Pauli matrices is denoted 
by cr, and we have absorbed a factor of 1/2 into the cou- 
pling constant Jk- Note that we have included a nearest- 
neighbor interaction in H Hubbard to qualitatively account 
for the poor screening of the Coulomb interaction from 
the insulating environment. The Coulomb-driven on-site 
and nearest-neighbor coupling constants are typically re- 
pulsive U,V > 0. However, in what follows we will treat 
these parameters as effective (phenomenological) ones, 
and assume that they include all possible contributions 
and renormalizations coming from the interaction be- 
tween stripe electrons and the nonmagnetic degrees of 
freedom of the environment, such as the electron- phonon 
coupling, or coupling to other electronic subsystems in 
the environment. As implicit in Eq. Q), we use the 
convention that the transverse coordinates take values 

= = 1 on the A- and B-arrays adjacent to the 
stripe. When convenient, we use the compact notation 

= S r j { i) (i = A, B) for the spins on these arrays. By 
assumption (vi) in Sec. I, the antiferromagnetic Kondo 
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lattice coupling Jk is weak, i.e. Jk "C Jn,t, allowing us 
to keep only the low-energy sectors of the stripe and the 
antiferromagnetic domains when analyzing its effect. 

The model in Eq. (JTJ is a modified version of that in 
Ref. 16 by having a coupling of two semi-infinite 2D an- 
tiferromagnetic domains to the stripe, one on each side 
of it, in this respect mimicking the geometry "seen" by a 
real stripe. As a consequence, with the assumption that 
the Neel order of the A and B domains are it— phase 
shifted relative to each other, we will be able to treat the 
metallic as well as the Mott insulating case (half-filled 
Hubbard band) within the same formalism. This is dif- 
ferent from the model in Ref. 16, where the assumption 
of a metallic stripe was crucial. Note that by turning 
off U and V in and keeping only one array of lo- 
calized spins, say, in the yl-domain, the Hamiltonian in 
([T]) collapses to the one-dimensional Heisenberg-Kondo- 
lattice model (HKL) 51 . This model has recently attracted 
a great deal of attention 2 ^2^ and we shall connect back to 
it when discussing our results in Sec. IV. Wc should here 
emphasize that having " built in" the presence of stripes 
into the model, we cannot address the issue of what actu- 
ally triggers the stripe formation. For this, one must turn 
to other approaches, some of them mentioned abov ofiffi 11 . 
A recent attack on the problem of stripe formation has 
made use of a spin-fermion model — which has the same 
Hamiltonian structure as l(T]l. but with the difference that 
there is no constraint on the mobility of the electrons in 
(J2J). In other words, the doped holes are now free to hop 
around on the 2D lattice. Monte Carlo simulations on 
the model suggest that the holes self-organize into one- 
dimensional charge stripes separating insulating spin do- 
mains (phase shifted by it across the stripes). This is the 
starting point when writing down our model in Q). 

Given the Hamiltonian in ljT)l. its partition function 
can be written as a path integral, 

Z = J V[c]V[^\D[n A \D[n B ]e- s[c ^ c ^ A ^ B] , (5) 

with a Euclidean action S. The electron operators are 
here simulated by Grassmann numbers (cj. a) c rta ), while 
the role of the localized spin operators S r ,■(.) are played 
by vectors Sfl r M) (i = A, B) which parameterize states 
in a coherent spin representation^,. Note that we have 
here used the short-hand notation c = {c r . a } and fti = 
{il r j(i)},i=A,B. In the limit of large spin, S — > oo, only 
diagonal matrix elements of a spin Hamiltonian survive 
in this representation. This makes the large-S" coherent 
state representation an efficient tool to mold a quantum 
partition function into a path integral. To retain quan- 
tum effects, however, present for physical values of the 
spin (S = 1/2 in the case of the cuprates), it is crucial to 
keep also non-diagonal matrix elements in the construc- 
tion of S. These produce a sum over Berry phases and 
contain a memory of the intrinsic quantum nature of the 
spins. The procedure is standaro-s, and one obtains the 



action 

s = J f \ T (^2cl a d T c r , a + h(c\c, sn Al sns^j 

+ iS J2 J2 (6) 

i=A,B r< jW 

where r corresponds to inverse temperature so that < 
t < (3 and the spin (Grassmann) fields are periodic (an- 
tiperiodic) in r. For the purpose of studying the low- 
energy dynamics we confine our attention to the zero- 
temperature limit (/? — > oo). The third term in (JBJ is 
precisely the sum over Berry phases 

$ rs = & dn r . ■ A(n ri ), (7) 
Jr r . 

one for each spin attached at site (r,j^) = r,-, i = A, B. 
Here Ttv is the closed loop traced out by $lr % in the in- 
terval [0,(3 — > oo], with A(f2 r J = (1 -cos flj) (sin 6>i)~V'i 
at each site r% playing the role of a vector potential of 
a unit magnetic monopole located at the center of the 
sphere = 1, parameterized by the spherical an- 

gles Oi and 4>i. The "instantaneous" Hamiltonian term 
i?(ct, c, SSI a, SSIb) in @ acts at (imaginary) time slice 
t and is obtained from by replacing electron and 
spin operators by the corresponding Grassmann fields 
(cl a ,c r7a ) and classical vectors (Sflr A , Stt rB ), respec- 
tively. 

Next, to obtain the low-energy continuum version of 
the Hamiltonian in we shall first review the standard 
constructions for a Hubbard chain and a 2D Heisenberg 
model, and then elaborate on the more intricate Kondo- 
lattice interaction which couples the two subsystems. 

1. ID Electron Chain: Low-Energy Theory 

The low-energy (field theory) approach to the ID (ex- 
tended) Hubbard model in © is based on the assumption 
of weak electron-electron interactions. Thus, assuming 
that \U\ , \V\ <C t, we linearize the spectrum around the 
two Fermi points ±kp (hp = n e Tr/2ao, where oq is the 
lattice spacing and n e is the electron density), and de- 
compose the original lattice operators into right-moving 
(i? Q ) and left-moving (L a ) chiral components: 

Cra e lkFrao R a (r) +e- lkFra «L a (r) 

- v 7 ^ [e lkFX R a (x) + e- ikpX L a (x)] , (8) 

where in the second line we have taken the continuum 
limit rao — > x. Defining local charge and spin densities 

Jr = 'R a Ra'> Jl ='L' a L a :, (9) 

Jr = ■ ■^RaO'apRp', J l ='--L) x cr a pLp:, (10) 

with repeated spin indices summed over, and with the 
normal ordering : ... : taken w.r.t. the ground state of 
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the free system, it is now straightforward to write down 
the low-energy continuum version of the U — V Hub- 
bard chain in (J2J. The weak interaction preserves the 



important property of spin-charge separation, and one 
can write the theory on the form H Hubbard = H c + H s , 
whereS^: 



H c = J dx{:J R J R : + :J l Jl- - 9q c Jl Jr - 9ouhn e {R\R\L-[Li +h.c.)}, 



(11) 



H s = 2ttv, 



lx{ : 



jz jz . 



■JiJi- 



■ 9osJIJr 



9o±(J x L JR + J y L J y R 



>}■ 



(12) 



r 



The velocities of the charge (c) and spin (s) excitations, 
governed by H c and H s respectively, are given by 



Vf ■ 



a {U + 6V) 
2tt 



vf ^ , (13) 



with vf — 2aotein(7rn e /2) the Fermi velocity. The small 
dimensionless coupling constants in and l|12|) are 
given by 



9oc = -a (U + 6V)/ttvf, 
9ou = -a (U - 2V)/itvf, 
9as = 9o± = aa(U - 2V)/irv F - 



(14) 



The Kronecker delta multiplying the Umklapp term in 
(|llfl signifies that this term survives the phase fluctua- 
tions (originating from the chiral decomposition in (JSJ) 
only for a half-filled electron band (n e = The trans- 
verse component of the spin current coupling in lJT2*)l. 



Here 7r ; = ±1 is the parity of the sublattice to which 
the site r t belongs, the unit vector n (suppressing the 
coordinate n for ease of notation) represents the local 
direction of the Neel order parameter field, and ao£ is 
a small orthogonal ferromagnetic fluctuation component, 
i.e. \a £\ < 1, with n • £ = 0. The phase factor (~1) 5 ' B 
in 11611 appears because from now on we take the Neel 
fields in the A and B domains to be 7r-shifted relative to 
each other. (The choice of reference vector in the stag- 
gering factor can be made arbitrarily, with no effect on 
the physics.) We here note that for a site-centered stripe 
embedded in a spin-1/2 environment current estimates 
predict that this is a viable assumption for stripe elec- 
tron densities n e < 0.6 in the limit where Jk ~ Jj^S*S. 
However, for Jk <C Jh, as assumed here, the critical 
density is expected to be larger. 



9o±{JlJr 



--^(Rl^Ll^+h.c), (15) 



describe backscattering of electrons. We have marked the 
spin-sector parameters in <|12fl , (|13|l and (|14|l by a " tilde" 
as a reminder that these will be modified when coupling 
the stripe to the localized spins in the environment. How 
this comes about is discussed next. 



2. Phase - Antiphase Antiferromagnetic Domains: 
Low-Energy Theory 

In the presence of antiferromagnetic correlations in the 
insulating domains - as seen experimentally in the stripe 
materialsi - the partition function in (J3J) at low energies 
is dominated by paths with 



i 



a$&(ri)n{ri)+aol{ri). (16) 



With n a slowly varying smooth field, (|16l) spells out 
the assumption of finite-range antiferromagnetic order. 
We should stress that n and i are taken to be indepen- 
dent fields, constrained only by the orthogonality condi- 
tion. This implies a doubling of degrees of freedom, which 
in principle should be corrected for when regularizing the 
theory. However, for our present purpose, to pinpoint the 
leading instabilities of the stripe electron dynamics due to 
the interaction with the antiferromagnetic domains, this 
issue is immaterial. The normalization of $7 in i|16|) is 
only preserved up to O(oq), but this is sufficient since we 
are interested in the long- wave length limit. In this limit 
we let the lattice spacing ao in the a;— direction (parallel 
to the stripe) go to zero, expand all terms in the action 
S which contain the spin fields up to 0(ag), and then do 
the replacements qp y^ r — > J dx, n ri {T) —* vim) (t, x), 
and &ri{i~) — > ijtt) (t, x). The result is a field theory for 
the independent orthogonal fields n and £, with an action 
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Si[n,£\ = - / g ^ I dx ( dr — {n j{i)+1 {T, x) - n j(l) (r, 
_,(*) 



CO ( T > 



5^ y"c?x J dT(n j(l) (T,x) x d T n j(l) (T,x)^ ■ £ j(t )(T 7 x) + S hB erry[n], i = A,B, 



(17) 



with 

poo 

Si,Berry[n\=iS^2"fr z (-l) S ' B / dr A(t, ri)-d T n(T, r 2 ) . 
r Jo 

(18) 

In the standard approach to the 2D antiferromagnet 53 
one would now integrate out the rapidly fluctuating £- 
field from l|17f) . After taking a continuum limit in the 
y— direction this would produce the familiar nonlinear a 
model (NLerM), with an added sum over Berry phases, 
describing the slow long- wavelength dynamics of the Neel 
order-parameter field n. In the present case, however, 
the localized spins adjacent to the stripe enter also in 
the Kondo lattice interaction and this must be taken 
into account before one attempts to integrate out the 
£— field. This problem is addressed next. 

3. Kondo Lattice Interaction: Mean-Field Decoupling 

Using © to write the electron spin density in Q 
in terms of the chiral fields, and then replacing the 
spin operators S r jd) in by the corresponding vec- 
tors Sfl r _j(i), decomposing these as in (II 61) . we obtain, 
expanding to 0{oq), 

Hliondo = + H n , (19) 

where 

H £ = J K Sal Ar • (£ ( r A) + £ { r B) ), (20) 



H n = J K Sa ^(-l) r A r ■ (n^ + n^), (21) 

r 

with 

A r = [e 2lkFra °Ll a R riP + e- 2lkFrao Rl a L r ^ 

+ Ll >a L TtP + Rl <a Rr >0 ]a al3 , (22) 

measuring the spin density on the stripe. We have here 
used the notation introduced after Eq. (0} , implying that 

£^ = £ r j(i) = i and = n r j(») = i. By the assumption 
that the Neel order directions of the A and B domains 



are shifted by 7r relative to each other it follows that 
ni^ = — rtr B \ and thus Hn vanishes. It is here im- 
portant to realize that there is no relative staggering of 
the ferromagnetic A and B components in J3UJ). These 
rapidly fluctuating fields are independent, with no corre- 
lations across the stripe. This leaves us with Hg, which 
in the continuum limit, using (|10|) with J = Jr, + Jr, 
takes the form 

F^M,*./«(.).(^(,) + < W W ). (23) 

We have here dropped the nonchiral terms that mix L- 
and i?-fields since these are washed out by the rapid 
phase oscillations in l|22l) . 

The fact that the stripe electrons couple manifestly 
only to the fast ferromagnetic components of the local- 
ized spins — independent of whether the stripe is metal- 
lic (Hubbard band away from half-filling) or insulating 
(half- filled band) 58 — does not mean that the Neel order 
parameter dynamics is completely decoupled from the 
stripe electrons. The Neel field re-enters the problem via 
the orthogonality condition n£ — which constrains the 
ferromagnetic component to a plane that follows its slow 
and smooth fluctuations. As we shall see next, part of the 
interaction in l|23|) can be absorbed as an effective spin 
density interaction among the stripe electrons. Since at 
low energies the Neel order direction is essentially con- 
stant over large patches in Euclidean space-time, this in- 
teraction will effectively be pinned in spin space, and 
hence break the spin rotational invariance of the electron 
spin dynamics on the relevant time- and length scales. 
This symmetry-breaking effect, driven by the Neel or- 
der in the environment, will dramatically influence the 
correlations of the stripe electrons. 

For simplicity, we now treat the A and B domains 
separately. Starting with A, and collecting all terms in 
the action containing the field defined on the — 
1 spin array that couples to the stripe, we find from i|17|) 
and H23f) the contribution to the partition function: 

= J V[£^] c- 5 ^'] , (24) 

with 
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S[£ (A} } = J dr Jdx 



4J H aoS 2 (£ (A) ) 2 + 2J K Sa (J L + Jr) • £ (A) - iS(n^ x d T n^) ■ £ {A) 

I 



(25) 



Note that we have again changed to the compact notation 

, introduced in Eqs. I|20|l 
and H21|) . The integral in (|24|) is Gaussian and can easily 
be carried out. We obtain 



exp [ - / drdx 



(wfr- 1 ^ =exp(-5'^), (26) 



where T = 4J H a S* 2 l, w = (2J K Sa )J± - iS^n^) x 
<9 r rc/ A )). We have here defined J± = J — (J-n^)n^, 
with «7 = «7l + «7/j, as the piece of the electron spin den- 
sity that — via the constraint ■ £^ = — survives 
the projection onto £^ A \ Thus, from (j26|) we obtain the 
effective action coming from fluctuations in 



ceff 



dr / dx 



a o J K r 



J 



iJ, 



K 



8irJ H 



(Ah 



1 



WJhcio 



(27) 



r 



Let us in turn discuss the different contributions to ^f/j : 

First term in \2l\j . The first term is an anisotropic spin 
interaction among the stripe electrons, induced by the 
coupling to the Neel-ordered spins in the A— domain. 
As we have already noted, this interaction follows the 
slow fluctuations of the n( A )-field along the stripe, with 
Jx constrained to a plane orthogonal to . To make 
progress we shall treat nS A ' in a mean-field formulation, 
and take it to be in a fixed (but arbitrary) direction (n) , 
defined by the antiferromagnetic order in the ^-domain. 
Introducing a coordinate system (x,y,z) with z in 
the direction of (n), and using the operator identity 
Jl/bJl/r = \ J L / R ' J l/Ri valid for chiral bilinears, the 
first term in 127|) is then seen to add the interaction 



Hint — 



2Jh 



J dx (: JiJi : + : J H J H : +J£JJj + JfJ H ) 

(28) 

to the spin-sector stripe Hamiltonian in (|12|l . The terms 
diagonal in L and R in H28|l are forward scattering 
terms which renormalize the effective spin velocity v s on 

the stripe, v s — > vi A ^ = v s — ^^/(AttJh), while the 
terms mixing L- and i?-fields describe backscattering 
of electrons and hence, when added to i|12[l . shifts the 
corresponding coupling: go± — > g± + aoJx/(ATVVs Jr). 

Second term in \2 r l\) . Let us first recall*? 4 that 
ii( A ) x d T n,( A ) is an angular momentum density, which, 
at the extremum of the action in l|17[) is locked to the 
ferromagnetic component: x <9 T «/ A ) ~ £^ A K Since 

is a rapidly fluctuating field, the second term in 
(|T7Jl . being a pure phase, for this case averages to zero 



already on finite patches in Euclidean space-time, and 
will be ignored in the low-energy limit considered here. 
This amounts to neglect fluctuations away from the 
cluster of paths that dominate the action (|17(l for the 
localized spins when decoupled from the stripe. As 
Jk <C Jh, we do not expect these paths to change much 
when inserting the stripe, and the argument applies 
also in the presence of the stripe. We shall discuss the 
limitations of this mean-field type argument below. 

Third term in \27j) . The last term in (|27|l . con- 
taining only the Neel field and its time derivative, 
should be assembled with the spin action in (|17|l . Then, 
integrating out all£^A)- fields from ||SJ - in exact analogy 
with the one-dimensional treatment of £j(A) = i in (|26|l - 
taking a continuum limit in the y-direction, and using 
(n x d T n) 2 = (d T n) 2 , one obtains the effective action 
for the order parameter field in the A-domain: 



S[n] 



I 

2 5o Jo 

Sphase [I. 



dr 



dx 
x Jo 



dy 



c Vn 2 +- — 

C OT 



(29) 



where the first term is the action for a nonlinear a 
model (NLcrM) in a semi-infinite plane, with parameters 
Qq 1 = S/V8a , c = ^f&JSa^. One piece of the original 
sum over Berry phases in (jf 7|) has been absorbed in 
the NLcrM, while the part H18|) containing the Neel 
field only, is left as a global phase S p h aS e[n] in (23 - 
This phase is an alternating sum over the solid angles 
$[n(ri,r)] swept by the local n(ri, r)-fields as r goes 
from to co. As long as there are no disordering or 
finite-size effects causing discontinuities in the Neel field, 
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Sphase [n] will be averaged outiS. For this reason we will 
ignore it for the moment. For the more realistic case 
of a finite-width antiferromagnetic domain, modeled, 
say, by a spin ladder—, S p hase[n} will come into play, 
requiring a more careful analysis. We shall return to 
this important issue in Sec. III.D. 



The analysis carried out for the A-domain above can 
be repeated step by step for the S-domain, and the fluc- 



J 



tuations in € ( - B - ) are seen to give a contribution identical 
to that in |J27J|. with the index "A" replaced by "B". 
Summing the contributions from the two domains, it fol- 
lows that the stripe electrons get described by an effective 
low-energy Hamiltonian 



B stripe B c -\- B s 

with H c defined in Eq. (H) , while 



(30) 



H s = 2ttv 



J dx{ : Jt Jt : + : J* R J' R : -g 0s JU R - g^{J x L J R + J*J&}, 



(31) 



r 



where 



gos 



ttvf 



(U-2V), 



ttvf 

Vs 2vr J H ' 



T 2 

,J K - 

2J H ' 



(32) 



with v s defined in Eq. (|13() . 

Thus, in contrast to the charge degrees of freedom 
which remain untouched by the coupling to the antiferro- 
magnetic environment, the spin dynamics on the stripe is 
strongly renormalized by this same interaction and gets 
controlled by an effective U (I) -symmetric spin Hamilto- 
nian H31fl . As for the low-energy processes in the decou- 
pled antiferromagnetic environment, these are described 
by two independent NLc models, one for each domain 
i = A,B, as defined in H29|) . (It may be worthwhile 
pointing out that as long as the Neel field is protected by 
the low-energy thermodynamic limit, the phase Sphase[ n ] 
in H29|) remains inactive^.) 



B. Diagonal Stripes 

The construction of the low-energy theory for a stripe 
running along the diagonal of the square lattice (FIG. 2) 
closely parallels that for a collinear stripe in the previous 
section. 




FIG. 2: Diagonal stripe structure. 



We again model the isolated stripe by an extended 
U — V Hubbard chain, but with a new hopping matrix 
element t' < t (since the overlaps between tight-binding 
orbitals along the diagonal of a lattice plaquette is ex- 
pected to be smaller than along the bonds). As a result, 
the continuum theory in and l|12l) still applies, but 
with t replaced by t' (implying a shift of the effective 
velocities defined in 11311 1. As for the antiferromagnetic 
domains, we expect that the order parameter dynamics is 
still described by the NLctM (|29|) in the bulk (away from 
the stripe): By assumption, the coupling to the stripe is 
weak, Jk <C Jh , and can only perturb spins in its imme- 
diate neighborhood. Here we are interested in the reverse 
effect, and will explore how the Kondo lattice interaction 
(0J affects the electron dynamics on the diagonal stripe. 

For this purpose, let us isolate one array of localized 
spins adjacent to the stripe, say the A-domain. We label 
the corresponding spin operators S^ A \ where r is a lat- 
tice coordinate running along the stripe that labels also 
the horizontal lattice axes that pierce the stripe at the 
corresponding sites. 

The Sl AS> spins interact with their neighboring spins 
on the parallel array, call them, via the terms 

Barray = Jh ^ ] (*^r ^ ' ^r— 1 ' )• (^3) 

3 



It is clear from the geometry that this interaction will 
generate an effective ferromagnetic coupling between 
the nearest- neighbor -spins, induced by a double- 
exchange via the two neighboring S^-spins. As seen 
in FIG. 2, this means that the local ID magnetic envi- 
ronment sampled by the stripe electrons via the Kondo 
exchange is ferromagnetically ordered. Does this imply a 
different induced interaction among the stripe electrons 
as compared to the collinear case in (|31fl ? To find out, 
let us first write down the full Kondo lattice interaction 
for the diagonal stripe, including the neighboring array 
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of S^, ^-spins from the B-domain: 

HKondo 

= JkJ2 c L<7«pc rJj ■ (S^ + sii\ + S< B > + S { r B \), 

r, a 

J K >0. (34) 

Following the same route as in Sec II. A. 3, using the de- 
compositions (JBJ and (|16J) . we obtain: 

H Rondo = Hg + H n , (35) 

where 

H t = J K Sal Ar • (£ { r A) + 4+i + ^ + (36) 

r 

H n = JkSclo A r ■ («£ A) +«r+i +^ B) +«r-i)i (37) 

T" 

with A r the local spin density operator on the stripe 
defined in (|22(l . Similar to the collinear stripe treated 
above, the 7r-phase shifted Neel order across the diagonal 
stripe - along horizontal as well as vertical directions - 
implies that ~ — ni B \ raj^H = —n^\ (see FIG. 2), 
and it follows that Hn in 13 7|) vanishes, independent of 
the value of the stripe band filling n e /3S. (Note that, in 



contrast to the collinear case in i|21[) . there is no stagger- 
ing factor in Eq. H37J1 : the local magnetic environment as 
seen from the stripe is uniform along the stripe.) We are 
thus left with in 13(i|) . Taking a continuum limit 

a ^^ dxdyS(x-y), (38) 

r 

and doing a gradient expansion to O(ag), we obtain 

4J K Sa J dxJ{x,x) ■ (e {A) (x,x)+£ (B) (x,x) S j , 

J(x,x) = J L (r)+J R (r), (39) 

with the "diagonal" stripe coordinate r replacing (the 
implicitly defined variable) x in the definition of Jl/r 
after Eq. I|23|l . Given l|39|l . we now again focus on the A- 
domain, and isolate the piece of (|39|l containing only the 
£^-field, Hg(A) call it. Inserting a time-dependence and 
letting Hg(A) act at (imaginary) time slice r, integrating 
over the slices, and adding the result to the ordinary 
semi-classical spin action Sa [n,£] for the A-domain (cf. 
with Eq. © for the corresponding collinear case), we 
obtain 



/■oo 

S A [n,£] + / drH £ 
Jo 





OO POQ pOO 

dx dy dr 

oo J x JO 



JhS 



^- ((Vn^)) 2 + 8(£^) 2 - iS(nW x d T n^) ■ *< A >) + 4J K Sa Q J ■ £ {A) S(x - y) 



(40) 



Writing down the full A-domain action in l|40l) we have 
ignored the presence of the sum over Berry phases (cf. 
again Eq. © for the corresponding collinear case), since 
it does not involve the £-field and hence does not couple 
directly to the stripe electrons. Also note that compared 
to our treatment of the collinear case in Sec. II. A. 3 we 



have here shortcut the analysis by taking a continuum 
limit in the y-direction before integrating out the £-ficld in 
Q40JI. Carrying out the integration we obtain an effective 
A-domain action <S^f/) generated by fluctuations in the 
ferromagnetic components of the localized spins: 



qeff 



dx 



dy 



dr 



0-0 J K 

Jh 



iJk 

Jx • Jx + ^J~( n x dm) ■ Jx ) 6(x-y) + 



(n x d T n) 



(41) 



with Jx defined after Eq. I|26|) . Doing a saddle-point in (|41|l . we obtain — in exact analogy with the collinear 
approximation and dropping the rapidly oscillating phase 
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case — the induced stripe-electron interaction 



Hint — 



2aoJ g 
Jh 



dr(J x L J x R 



(42) 

with r the "diagonal" coordinate along the stripe. 

This Hamiltonian is almost an exact copy of that for 
a collinear stripe, cf. Eq. I)28[l. with the only difference 
that the magnitude of the coupling is larger by a factor of 
4. (This can be traced back to the fact that on the lattice 
a diagonal stripe electron couples simultaneously to two 
localized spins in the A-domain.) Adding the identical 
contribution from the B-domain, it follows that the low- 
energy spin dynamics on a diagonal stripe is described by 
the same effective Hamiltonian <|31[) as for the collinear 
case, but with the renormalized parameters in l|32|) mod- 
ified by Jk — ► 2Jft-- (As mentioned above, the hopping 
matrix element t' is also different from that which enters 
the free part of the Hamiltonian for a collinear stripe; cf. 
Eq. (|3U[1 with vp — 2aofain(7m e /2). However, as this has 
no impact on the problem studied here, from now on we 
use a single label t for both types of stripes.) 

In contrast to our analysis for the collinear stripe con- 
figuration, we have here not attempted to derive the full 
effective spin action for the decoupled environment. As 
we have already noted, away from the stripe the Neel or- 
der parameter dynamics is described by a NLcrM with an 
added Berry phase, as in (|29|l . By inspection one finds 
that close to the diagonal stripe the Berry phase gets in- 
fluenced by the unusual boundary condition associated 
with the diagonal stripe orientation. Thus, our results 
- here derived for a semi-infinite 2D geometry - may be 
of limited applicability for the case of "diagonal" finite- 
width or "spin ladder" environments (see Sec. III.D). 
Their study is an interesting problem, but we here leave 
it for the future. 

The fact that the same effective interaction appears 
for the diagonal and collinear stripe structures (up to 
the trivial Jk <-> 2Jk shift) reflects its origin in the 
coupling of the electronic spin density to the uniform 
t- components of the localized spins. These components 
are confined to a plane orthogonal to the Neel-direction, 
and are blind to whether the local Neel-field adjacent 
to a stripe is staggered (as for a collinear stripe) or 
uniform (diagonal stripe)Si. The coupling constant 
J K I Jh embodies the second-order process that drives 
the induced interaction between the stripe electrons: An 
electron exchanges spin with the environment (~ Jk) 
and another electron arriving at the same lattice site flips 
back the localized spin by a second exchange (~ Jk), 
resulting in an effective spin exchange (~ J K ) between 
the two electrons. Since only that part of the spin 
exchange that couples to the ^-components survives, the 
effective interaction becomes anisotropic. The 1/Jh~ 
dependence of the process is also expected: The larger 
the spin stiffness of the antiferromagnetic environment, 
the smaller the probability for the double exchange to 



C. Effective Stripe Hamiltonian: Bosonization 

As we have seen, the low-energy electron dynamics on 
collinear as well as diagonal stripes - taking into account 
a weak spin exchange with the environment - is described 
by an effective Hamiltonian 130f) . with the amplitudes for 
forward- and backward scattering renormalized by the 
exchange. This Hamiltonian embodies a spin- anisotropic 
interaction among the electrons, well-defined on length- 
and time scales over which the environment is magneti- 
cally ordered. To analyze the consequences for the stripe 
electron dynamics we shall use Abelian bosonization to 
map the model onto two independent quantum sine- 
Gordon models (in the weak-coupling limit)- one describ- 
ing the collective charge excitations, the other the spin 
excitations - and then perform a renormalization group 
analysis to identify the leading instabilities of the system. 

The method of bosonization is well reviewed in the 
literature£2i£ii§£, and we here only sketch the most im- 
portant steps so as to fix notation and conventions. The 
standard bosonization formulas for spinful chiral elec- 
trons are given by^S: 



R a (x) 

L a (x) 
40*) 



1 



%/27ra 
1 



n a e 



y/^TTOQ 



1 



n a e 



(43) 



V27TO0 

1 



V27TO0 



n a e 



n a e 



-iv4Tr<f>L,a(x) 



(44) 



Here 4>R, a (x) and <f>L.a{x) are right- and left moving 
bosonic fields respectively carrying spin a = ],[. The 
Klein factors r\ a and fj a are inserted to ensure that the 
anticommutation relations for electron fields with differ- 
ent spin come out right. They are Hermitian and satisfy 
a Clifford algebra 

{na^p} = {n<x,m} = 2 <Wh {v<x,vp} = Q- (45) 

One next introduces charge (c) and spin (s) fields (p CjS 
and their duals t? c , s : 



(^ + ^)/V2, 



#c = (ei + e L )/V2, (46) 

0. = (0 T - (47) 



where 



Then, using the identities 

Jr + Jl = 
Jr — Jl = 



-(l/y/n)d x ip 



C ! 
C ) 



(48) 



(49) 
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Jr + J z l 



T z — 
J L — 



-(1/V^7r)d x cp s , 



(50) 



together with i|43|) and (|44|l . we can translate H str ipe 
H c + H s in l|30|) into bosonized form: 



dx 



t (Xr\ J 



(51) 



2m s 



-K COS 



(V8^7^.)}. 



(52) 



We have here introduced the rescaled charge and spin 
fields 



0' 



Kl'fa,., (53) 



and the short-hand k = VlVlV^Vi- To leading order in 
the coupling constants the sine-Gordon model parame- 
ters -Kc(s) and m c / s -\ are given by 



2{K C - 1) 



,9oc 



- — (U + 6V), 
ttvf 



27rm c = -g 0u = -(U-2V), 

TTVp 



(54) 



in which the above product k — 77^- 77^ 77-p r)i of Klein factors 
has the form 



K = 1 



1 (J 3 . 



(60) 



This matrix is diagonal with eigenvalues ±1. Provided 
that all relevant correlation functions to be calculated 
contain only products of Klein factors which are simulta- 
neously diagonal with k we can pick the eigenstate with 
eigenvalue +1, say, and then ignore the rest of the Klein 
Hilbert space (allowing us to do the replacement k — ► 1 
in """J} and (152*)) ). We will come back to this point below. 



2(K s -l) = 90s = -^(U-2V), 
ttvf 



27rm, 



9ox=— (U-2V + (3^), (55) 



with 



Vc = v F + ^-(U + 6V), 
a 



T 2 



v F -^(U-2V + 2^), 



where 



= 



f 1/2 collinear stripe 



2 diagonal stripe 



(56) 
(57) 

(58) 



Note that in obtaining (|51[) and (|52(1 . terms corre- 
sponding to scattering processes which lead to a renor- 
malization of the Fermi velocities in second order in 
the couplings, as well as strongly irrelevant terms ~ 
cos(y / 8n K c (p' c ) cos(y/8Tr K s ip' s ) describing Umklapp pro- 
cesses with parallel spins, have been omitted. 

The product of Klein factors in (|51() and (|52(l acts on 
a Hilbert space different from the boson Hilbert space 
and introduces a certain ambiguity into the formalism. 
We resolve it by choosing a representation of the Clifford 
algebra in terms of tensor products of Pauli matrices and 
the identity operator, 



'7r 



01 
02 



cti , 



<7 3 (g> Cl , 
1® <7 2 , 



(59) 



III. PAIRING AND DENSITY CORRELATIONS 

A. Renormalization Group Analysis 

The mapping of the effective stripe Hamiltonian in 
Eq. I|30l) onto the quantum theory of two independent 
charge and spin Bose fields, manifestly shows that the 
collective low-energy charge- and spin dynamics on the 
stripe remains separated in the presence of a magnetic 
environment. This allows us to extract the ground state 
properties of the stripe electrons by performing indepen- 
dent renormalization group (RG) analyses of the charge- 
and spin sector sine-Gordon Hamiltonians. The RG flows 
are of Kosterlitz-Thouless type, with effective coupling 
constants gi (i — c, s, u, _L), governed by the equations 6 - 3 



dg c /dl = -g 2 u , 
dg u /d£ = -g c g% 



for the charge sector, and 



dg s /dt = -g 2 ± , 
dg±/d£ = -g s g±, 



(61) 



(62) 



for the spin sector. Here £ = £n(a/ao) with a a renormal- 
ized length, while gi (£ = 0) = goi are the bare parameters 
that enter Eqs. ("5"4"}) and J55jl. We shall denote by K c ( s ) 
and m c ( s ) the corresponding renormalized sine-Gordon 
parameters connected to gi via the same Eqs. (|54H and 

(E""J. 
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The flow lines lie on the hyperbolas 

2 2 _ 2 2 fao\ 

9 c (s) ~ 9 U (±) — 9qc(s) - 9ou(±) ' \ 06 ) 

and — depending on the relation between goc(s) arL d 
3o«(±) ( or ? equivalently, the bare sine-Gordon parame- 
ters -Kc(s) and m c r s -\) — exhibit two types of behaviors 
(cf. FIG. 3): 




FIG. 3: Renormalization-group flow diagram for the spin 
(flsiflx) an( f charge (gc,gu) sectors. The arrows denote the 
direction of flow with increasing length scale. 

Weak coupling regime. 
When g 0c(s) > \g 0u (±)\ ( K c(s) - 1 > 7r|m c ( a )|) we are 
in the weak coupling (Luttinger liquid) regime: g u (i_) — » 
0, implying that the renormalized masses rh c (s) scale to 
zero. The low-energy, long-wavelength behavior of the 
gapless charge (spin) degrees of freedom is thus described 
by a free scalar field 

H c[s) = V ~^Jdx {(d x ^ c(s) f + {d x $' c(s) f} . (64) 

Ignoring logarithmic corrections^ coming from the slow 
renormalization of marginally irrelevant operators near 
the fixed line g u (±) — 0, the large-distance behaviors of 
the charge- and spin field correlators and their duals are 
given by 

(65) 

( e i \/2V^<( a )(^) e -» v ^7^<( S )(o)^ _ \ x yUK M . 

(66) 

Hence, the only parameters controlling the low-energy 
behavior in the gapless regimes are the fixed-point values 
K^fs) (Luttinger liquid parameters) of the renormalized 
coupling constants K c ^ »1| <? c (s)/2- 



Strong coupling regimes. 

When g 0c{s} < \g u(X)\ ( K c(s) - 1 < 7r|m c ( s )|) the 
system scales to strong coupling. The two sepa- 
ratrices <?oc(s) = ±|<?0u(_L)| divide the strong cou- 
pling regimes for charge and spin into two sectors 
respectively: (i) g 0c ( s ) < ~\gou(±)\, where the in- 
crease of ffc(a) an d |5«(-L)| i s immediate, and (ii) 
-|pOu(±)| < 9oc(s) < \gou(±)\> where one observes a 
crossover from a weak-coupling behavior at intermediate 
scales (<7c( s ) ~ |ff u (i)|) to strong coupling at larger scales 
(Sc(s) ~ — |g M (j_)|)^i. Depending on the sign of the bare 
mass m c r s ) in (|54|l and l|55|l . the renormalized mass 
fh c ( s ) is driven to ±oo, signaling a flow to one of the two 
strong coupling regimes, with a dynamical generation 
of a commensurability gap A c ( s ) in the charge (spin) 
excitation spectrum. The flow of |to c ( s J to large values 
indicates that the cos-term in the sine-Gordon model 
dominates the large-distance properties of the charge 
(spin) sector. With the cosine-term being the dominant 
one, the values of <p' c / s \ will tend to be pinned at the 
minima of the cosine potential. For m c ( s ) < these are 
at yf&nK c { a )<ff c r s \ = 27m, with n an arbitrary integer. 
Since <^(s) are an g u l ar variables one cannot distinguish 
between different n, however, and the "negative mass" 
condensation is defined by (tp' c r s \) = 0. Similarly, for 
m c ( s ) > the minima are at ^JSttK^^lp'^^ = irn and 

the fields order at ^Jtt /8K c ( s y To summarize, there are 
two strong-coupling regimes where the fields <p' c r s ) get 
ordered with the expectation values 

Wc(s); \0, m c(s) <0 • [0,) 

Note that the signs of the bare masses in (|6*7jl are 
contingent upon the choice of truncated Klein Hilbert 
space in Sec. II. C, where we have taken k —> 1 in l|51|) 
and l|52p. This has no effect on the physics, however, 
since a transposition of the two strong-coupling phases 
above (via the alternative choice n — * —1) would be 
followed by a subsequent redefinition of any relevant 
correlation function, thus producing the same value of 
any observable. 

Having exposed the properties of the weak- and strong 
coupling regimes, let us apply the results first to the 
charge sector. By inspection of the "bare" values of the 
coupling constants in the charge sector, Eq. 1)540. one 
easily finds, using Eqs. (|67|l . that for a half-filled band 
(n e — 1) this sector is gapped for U > 2\V\ and for 
U < 2V when V > (strong coupling regimes). In the 
former case rh c — > — oo, implying that 

(<p' e ) = 0, (68) 
while in the latter case m c — > oo, with 

{<p' c ) = V^Wc (69) 
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For any other values of U and V , but still at half-filling, 
we are in the weak- coupling regime, corresponding to a 
gapless charge excitation spectrum. The charge degrees 
of freedom are here governed by the free Bose field in 
Eq. (JB3J, with the fixed-point value of the charge param- 
eter 

iC-l + — VV(U + 2V)>1, n e = l, (70) 

TTV C 

obtained from with rh c = 0. The line U = 2V > 
corresponding to the fixed-point line m c — 0, K c — 1 < 0, 
is special. Here the low-energy properties of the gapless 
charge sector are described by the free massless Bosc field 
in l|t)4|) with K* = K c , reflecting its exact marginality. 

Away from half- filling (n e ^ 1) the bare mass term in 
(|5T1) is killed off for any values of U and V, and the charge 
degrees of freedom are described by the free Bose field in 
()64p. Analogous to the special line U — 2V > above, 
the correlations in H65J1 and Ijfifijl are now governed by 
the bare value of the Luttinger liquid charge parameter 
in 

K* = K c = l-^{U + 6V), n e + \. (71) 
Zirv c 

All of the above is familiar from conventional "g- 
ology" for Hubbard-type models 12 - 55 . Since the Kondo 
lattice interaction, Eq. J3J , does not couple to the charge 
sector the latter indeed behaves as if the electrons were 
isolated from the magnetic environment. 

Let us now look at the behavior of the spin sector, 
which is more interesting. As we have seen, the spin ex- 
change between the Neel-ordered environment and the 
stripe electrons breaks the SU(2) spin-rotational sym- 
metry in the effective theory. This implies that the 
spin sector is gapped for arbitrary U ^ 2V — PJ^/Jh 
(strong coupling regime). When U > 2V — (3J] ( /Jh 
the mass renormalization goes to +oo, whereas for U < 
2V — (3 Jft/ Jh the mass renormalizes to — oo. Reading off 
from 16711 , using (|55|l , this implies the spin field orderings 




Note that this result is independent of the band filling on 
the stripe. 

When Jk = and the stripe decouples from the en- 
vironment, the SU(2) invariance of the spin sector is re- 
covered, and the spin dynamics renormalizes along the 
separatrix g s = g±. As is well-known, the spin sector 
then gets controlled by the weak- coupling Luttinger liq- 
uid parameter K s ~ 1 + \g s — > K * = 1 when U > 2V 55 , 
whereas for U < 2V one stays in the strong coupling 
regime with (ip' s ) = 

Next, we want to exploit the RG results derived above 
to map out the groundstate phase diagram for the stripe 
electrons. We shall catalogue the different phases accord- 
ing to the values of (U, V, J\j Jn,n e ), and focus on the 



corresponding behaviors of density and superconducting 
pairing fluctuations. These are characterized by the cor- 
relations of the associated order parameters, which in the 
present case come in two guises: (1) conventional and (2) 
composite order parameters. Let us in turn review their 
definitions and bosonized representations. 

B. Order Parameters 

1. Conventional Order Parameters 

The conventional order parameters^iSiiSSiSl which 
may develop long-range correlations in this class of mod- 
els are those of short wavelength (2kpx) fluctuations 
of the site- and bond-located charge density, site- and 
bond-located spin density and superconducting singlet and 
triplet pairing. By using the chiral decomposition (|SJ) of 
the electron fields together with the bosonization dictio- 
nary in Eq. I|43|l - H5U|) it is straightforward to obtain 
their bosonized forms. 

• Charge density wave. 

This fluctuation is carried by the charge- 0, spin-0 exci- 
tations created by the operator 

Oc'Dw(r) = ^ 4.a c r,q 
a 

= ( Jr ^ ( r >^« (rWrLl (r)R a (r^i^ {r)L a (r)) 

(73) 

with S r — e 2lfcF ™°, and where Jl/r are the chiral charge 
currents defined in Eq. (0- Keeping only the finite- 
momentum modes (k — ±2kp) from the nonchiral terms, 
taking a continuum limit, and reading off from the dic- 
tionary (|43l) - (|50|l . one obtains the bosonized expression 

Ocdw(x) -> 

sin(V27r K c ip' c (x) - 2k F x) cos( v / 2TrK s ip' s (x)), (74) 

where we have used that 77-p 77-p and fj^rji are diagonal with 
the same eigenvalue on the truncated Klein Hilbert space 
chosen in Sec LLC 

• Spin-density wave. 

This is the simplest charge- 0, spin-1 vector order param- 
eter, and is defined by 

OsDw(r) = ^4,o° 'oipCrfi- (75) 

Bosonizing the x-component of 0(ji)sdw that creates 
k = ±2kp excitations, and dropping the trivial k = 
modes, one finds in the long-wave length continuum limit 

®Sdw( x ) ~* 

?7 T f7X cos(y/2wK c <p' c (x) - 2k F x) sin(V ' 2irK7 1 'd' s (x)). 

(76) 



15 



To obtain this form we have exploited the fact that the 
Klein factors fj^r]^, — fj^ , r^f/j and — r}]fj^ have the same 
action in the truncated Klein Hilbert space. (Note, how- 
ever, that 77-p 77^ is not diagonal on this space, and as a 
reminder of this we keep it explicitly in Eq. 176(1 ). In the 
same way one easily obtains 

fj^rji cos(-\/27r K c ip' c (x) — 2k F x) cos(v 2nKs 1 -d' s (x)), 

(77) 

and 

®Sdw( x ) ~~ * 

cos(y/2irK c tp' c (x) - 2k F x) sm(y/2wK s <p' s (x)), (78) 

where in the z-component we have put fj^Tj^ — ?7x?7x = 
constant. 

In the special case of a half-filled band (n e = 1) one can 
distinguish between the 2k F modulations of the charge 
and spin densities with extrema of the density profile 
located on sites or between sites - i.e. on bonds. There- 
fore at half-filling one should also consider order param- 
eters corresponding to the short wavelength fluctuations 
of bond-located charge and spin densities: 

• Bond-located charge density wave ( "dimer" ) . 



A dimerization instability is characterized by enhanced 
correlations among the charge-0, spin-0 excitations cre- 
ated by 

ObCDw(r) = y^c^qCr+i.c* + h.c). (79) 

ot 

Again, keeping only the k = ±2k F excitations, one ob- 
tains in the continuum limit 

ObCDw(x) — > 

cos{y/2wK c <p' c (x) - 2k F x) cos(y / 2irK s ip' s (x)). (80) 

• Bond-located spin- density wave. 

This is the vector order parameter that describes 
charge-0, spin-1 magnetic excitations centered on the lat- 
tice bonds: 

O bS Dw(r) = - y^(4.a cr afl c r+l,|3 + h.C.). (81) 
a.fi 

In the continuum limit the lattice shift in (|81|l shows up 
as an extra phase 7r/2 added to the ubiquitous phase 
2kpx (cf. the bosonized dimer operator in 1)801) ) and 
we thus identify the bosonized components of the finite- 
momentum part of ObSDw{x) as 



{rj-^rji sin (tJ2icK c <p' c (x) — 2k F x) sin [y/2^K s 1 i9' s (x)), i = x, 
fj^i s in (y/2T K^ip' c (x) - 2k F x) c os{\/ 2ttKs 1 $' s {x)) , i = y, (82) 
sin (y/2nK c ip' c (x) — 2k F x) sin (y/2irK s <p' a {x)') , i = z. 

Finally, we consider the two order parameters for (superconducting) pairing: 

• Singlet pairing. 

The charge-2e, spin-0 superconducting pairing modes on the stripe lattice are created by the operator 

o ss (r) = 4, T 4a = ^4w L l( r )+ <y r 1 4( r ) E l( r )+ L }( r ) iZ K r )+4( r )4( r )- ( 83 ) 

The k — ±2k F excitations produced by the chiral terms are the so called 77-pairing modest. The right-moving 77-pairs 
can be written as 

7] R (x) ee r\(x)B\(x) -> cxp{i\/2TrKc 1 v v c (x))exp(-i(y/2TrK c (p' c (x) - 2k F x)) (84) 

in the long-wave length limit, with the analogous expression for left-moving pairs, T) L (x) ee L^(x)L^(x). As these 
contain only the charge field and its dual, they are blind to the antiferromagnetic environment and hence we will not 
consider them here. This leaves us with the k — BCS singlet pairing operator contained in (|83|l . with the bosonized 
form 



Oss(x) mm exp(i\ / 2TrKc 1 K(x))cos(^2nK sV ' s (x)). 



(85) 
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• Triplet pairing. 

The charge-2e, spin-1 pairing modes are created by the lattice operator 

TS {r)=-icljcTCjy) a f j cl Js . (86) 

Again retaining only the k — modes, 

TS {r) -> -i 4 (r)( CT ^) a/3 4 (r) , (87) 
we obtain for the bosonized components in the long- wave length limit: 



exp^^Tri^V^x)) sm(y / 2irKi 1 ^ s (x)), i = x, 
O ts -> <j exp(z\/2^ c - 1 ^(cc)) cos(v / 2^ s - 1 ^(x)), t = 

ri^fjiexp(i\'2wK^'&' c (x)) sm(y/2TrK s <p' s (x)), i = z, 



2. Composite Order Parameters 

In addition to the conventional order parameters listed 
above we need to consider composite order parameters 
built from operators acting on the stripe electrons and 
the magnetic environment. The notion of composite 
order parameters was first exploited in the theory of 
superconductivity^*^, where it was realized that since 
any product of a particle-hole (i.e. charge-neutral) oper- 
ator and a Cooper pair operator possess charge 2e this 
composite can, in principle, describe some superconduct- 
ing state. By analogy, one may similarly construct com- 
posite CDW and SDW order parameters. 

We shall here introduce only composite order param- 
eters that may develop long-range correlations for the 
physically most interesting case of a stripe away from 
half-filling and with repulsive electron-electron interac- 
tions U — 2V > 0: composite CDW and composite singlet 
pairing. 

• Composite (site-located) charge density wave. 

A composite CDW order parameter is obtained by pro- 
jecting the conventional (site-centered) spin-1/2 SDW 
onto the difference between the localized spins on the 
neighboring A- and £?-arrays: 

O c ^cdw ~ Osdw ■ (S (A) - S (s) ). (89) 

Note that this expression is well-defined in the continuum 
limit for any stripe geometry: In particular, for the case 
of a diagonal structure, an electron at the r th site on the 
stripe couples to \ + in the A-domain, which 

in the continuum limit reduces to 2S^ A \x r ), dropping 
an irrelevant gradient term. Considering first a collinear 
structure, we need to keep only the staggered parts 
ni A \— l) r and ni B \— l) r of the localized spins since the 
correlations of the uniform components of Sl A ' and 
die out faslpS. With a phase-antiphase domain, as as- 
sumed here, we further have that = —n^ = n. 



Thus, from (|S3|l . the k = 2k F + 7r/a = 2k* F part of the 
composite charge density wave is given by 

Oc-cdw ~ Osdw • n(-l) r , (90) 

with r the stripe lattice coordinate, and with the 
bosonized components of Osdw written down in 
Eqs. (|76|l - J7HJ. It would be tempting to refer to the 
generalized Luttinger theorem^ to " explain" the appear- 
ance of the composite staggered CDW, Eq. (|90|l . As 
pointed out by Zachar— , the theorem asserts that the- 
ories belonging to the class of Kondo-Heisenberg lattice 
type models are expected to support a massless spin-0, 
charge-0 excitation of momentum k = 2k* F (reflecting the 
presence of a "large Fermi surface" due to the localized 
spins). However, in the present case the localized spins 
of the environment are assumed to be ordered (with the 
NLerM in Ij29(l describing the small fluctuations of the 
order parameter field n), and, as a consequence, time 
reversal symmetry - entering as a condition for the va- 
lidity of the theorem^ - is broken. Indeed, the case of 
a diagonal structure is different, and does not produce 
a k — 2k* F mode. Here the stripe electrons experience 
a local ferromagnetic environment (cf. FIG. 2), and the 
composite CDW now appears at k = 2k F (i-e. with no 
staggering) : 

0^c k D F w-OsDwn. (91) 

• Composite singlet pairing. 

By taking the product of the conventional triplet pairing 
operator Ots for the stripe with the difference of spin 
operators for localized spins, — S^, a composite 
singlet pairing operator can be formed as 

O c -ss - Ots ■ (S (A) - S {B) ). (92) 

For a collinear stripe this operator has two momentum 
components: a uniform k — composite singlet 
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with rapidly decaying correlations due to the incoherent 
fluctuations of and l^ B \ and a k = ir/ao staggered 
composite singlet 

o (fc=Va ) _ o TS -n{-lf (93) 

= -t J R+(r)( < T<r") a/ ,iJ(r).n(-l) r ) 

with r the discrete lattice coordinate along the stripe, 
and where n = n A = —n B . It is important to note that 

or 

™J a o) j g Q ^ unc jer time reversal (T : R <-> L, a — > 
— cr, re — » — re), as well as under parity (P : R <-> L), 
implying "odd-frequency odd-parity pairing" 29 . 

Turning to the case of a diagonal stripe structure the 
composite singlet pairing now occurs for k = (since 
the local magnetic environment appears uniform as seen 
from the stripe), and one has 

£> ( c-ss ~ -iRUr)(^) af3 Ll(r)-n. (94) 

Again, parity and time reversal are broken. It is impor- 
tant to point out that theoretical work«££ suggests that 
odd- frequency pairing is actually unstable for k = pairs 
(at least within Eliashberg-Migdal theory, where vertex 
corrections to the self-energy are neglected). This result 
becomes particularly intriguing when seen in the light of 
the diagonal — ► collinear stripe rotation associated with 
the superconducting transition observed in some of the 
cuprate a 38 i 39 (cf. our discussion in Sec. I): If the singlet 
pairing in the high-T c materials were of composite na- 
ture, the stripe rotation would precisely serve to stabilize 
the pairing by shifting the momentum from k = (di- 
agonal configuration with unstable pairing) to k = tt/ciq 
(collinear configuration with a stable, staggered compos- 
ite pairing mode)! 

In the presence of 2D Neel order, as assumed here, the 
re-field correlations are infinitely ranged in the ground- 
state, and the Cc-ss"" ano - ^c-ss operators may build 
up large-distance correlations that compete effectively 
with conventional triplet pairing. Whether this happens, 
and what other order parameter correlations may de- 
velop, will be studied next. For an extended discussion of 
composite order parameters for ID correlated electrons, 
we refer the reader to Refs. 26 and 28. 



C. Phases 

Equipped with the results in the two previous sections 
we shall now pinpoint the leading groundstate instabil- 
ities of the stripe electrons and list the corresponding 
phases: 



1. Half-filled band: n e = 1 

The phase diagram consists of five sectors: A, B, C, 
Dl, and D2 (see FIG. 4). 



• A phase: U > 2\V\. 

We include this case merely as an illustration of our 
formalism, as a half-filled band (one electron per site on 
the stripe) is somewhat special when combined with dom- 
inant repulsive on-site interaction. The reason is that one 
here expects to loose the phase- antiphase configuration as 
n e — ► 1, and instead recover the undoped antiferromag- 
netic state (with an in-phase Neel configuration across 
the "stripe"). In work by Zachar— , based on a stripe 
t-J model (corresponding to a "strong- coupling" limit 
Jk = Jh of our lattice model in (JTJ), it was suggested 
that there is a transition from the phase-antiphase to in- 
phase Neel configuration already at a band filling ~ 0.6 
(see also Ref. 57). Still, it is instructive to explore the 
assumption of a half-filled stripe embedded in a hypo- 
thetical phase-antiphase Neel background, and explore 
the consequences. At half-filling the charge excitation 
spectrum is gapped (A c ^ 0) when U > 2\V\. The stripe 
is thus insulating and the ordering of the charge boson 
with ground state expectation value (ip' c ) = suppresses 
the conventional CDW and superconducting correlations, 
but leaves behind the SDW and Peierls (dimerized) corre- 
lations. Turning to the spin sector, according to Eq. I|72|) 
there is a condensation at (ip' s ) = y/w/8K s . This kills off 
the Peierls correlations, and as the "in-plane" SDW 1 '"- 
correlations are seen to be incoherent we are left with 

{O z SDW {x)O z SDW {x')) ~ cos{2k F {x - x')) 

— > (-1)' x constant, (95) 

where in the last step we have reintroduced the discrete 
stripe coordinate x — ra , x' = (r + 0°o- Thus, given 
the hypothesis of a phase-antiphase spin background one 
obtains a long-ranged antiferromagnetic (Neel) phase for 
the stripe. The energy of this frustrated configuration 
grows linear with the length of the stripe and is hence 
unphysical, as anticipated. Note, however, that the ac- 
tual "in-phase" Neel configuration for n e = 1 does imply 
a long-ranged Neel phase for the stripe! 

• B phase: < U + PJ 2 K /J H < 2V. 

This phase is that of an insulator with a long-range 
ordered CDW: Both charge and spin excitations are 
gapped. The fields V'c(s) 01 "dered with ground state 

expectation values (ip' s ) = and (<p' c ) = ^tt/8K c respec- 
tively, and 

(Ocdw(x)Ocdw(x')) - (-1)' x constant, (96) 

with I defined after Eq. (|95|l . For the case of an isolated 
stripe (Jk = 0), results from weak-coupling perturbative 
renormalization group studies^S show that there is a 
continuous phase transition along the line U — 2V sepa- 
rating the SDW Z and CDW phases. The recent interest 
in the extended U — V Hubbard model was triggered 
by Nakamura^i, who found numerical evidence that for 
small to intermediate values of U and V, the SDW Z and 
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FIG. 4: The ground state phase diagram of the stripe electron 
system at half-filling. Solid lines separate different phases: 
SDW Z — long-range ordered spin density wave phase; CDW 
— long-range ordered charge density wave phase; SS — BCS 
(superconducting) singlet pairing phase; TS Z — triplet pairing 
phase (coexisting with a composite SS); bSDW z — long-range 
ordered bond-located spin density wave. As explained in the 
text, the ^-domain can only be realized with an "in- phase" 
Neel configuration across the stripe. The other stripe phases 
are assumed to coexist with antiphase Neel configurations of 
the localized spins. 



CDW phases are mediated by a bond-located charge- 
density-wave (bCDW) phase: The SDW 2 -CDW transi- 
tion splits into two separate transitions: (i) a Kosterlitz- 
Thoulcss spin gap transition from SDW Z to bCDW and 
(ii) a continuous transition from bCDW to CDW. An 
analogous sequence of phase transitions in the vicinity of 
the U = 2V line is an intrinsic feature of extended U — V 
Hubbard models with bond-charge coupling 66 . A similar 
effect is here caused by the Kondo coupling to the an- 
tiferromagnetic environment: At Jk ^ 0, along the line 
U = 2V only the charge gap closes. Therefore this line 
correspond to a metallic state with dominating antiferro- 
magnetic SDW Z and bSDW 2 correlations (since K c < 1), 
showing identical power-law decays at large distances: 

{ < ^ > SDw{ x )0 Z SDW {x)) = (®bSDw( x )®bSDw( x )) 

~ (-1)' x \x-x'\- K °. (97) 



• C phase: -f3J^/J H < U - 2V < and V > 0. 

Here again a charge gap opens. However, since in this sec- 
tor U — 2V < the bosonic charge field is now co ndensed 
with ground state expectation value (<p' c ) = y/n/8K c . 
This immediately leads to suppression of the site-located 
SDW Z correlations, and instead the bond-located bSDW z 
exhibits long-range order: 

{OtsDw{x)OtsDw(?')) - (-1)' x constant. (98) 



The line V = 0, U < is the crossover line from the 
insulating phases into the superconducting phases. On 
this line the charge sector is in the weak-coupling gapless 
(metallic) phase with K* = 1. However, the spin sector is 
massive along this line except a the point U = —(3J% / Jh, 
which marks the transition from a metallic phase at 
-PJ 2 K /J H < U < where the SDW Z , bSDW 2 and TS Z 
fluctuations show identical algebraic decay at large dis- 
tances 

(OsDwi^OsDwi 00 ')) ~ (^bSDwi^^bSDwi^)) 

~{O z TS {x)O z TS {x'))~\ x - x i\- 1 , (99) 

to a different metallic phase at U < — PJ^/Jh, where 
the SDW, bSDW and TS Z fluctuations are suppressed, 
while the conventional CDW, SS and Peierls correlations 
show identical large distance behavior: 

(Ocdw(x)Ocdw(x')) - (O ss (x)Oss(x')) 

)) ~ \x - x'\ . (100) 

This large degeneracy of metallic phases along the 
line V — is due to the SU(2)-charge ("pseudo-spin") 
symmetry of the half-filled Hubbard model 6 ^. This de- 
generacy is immediately lifted by an attractive nearest- 
neighbor coupling (V < 0), in support of superconduct- 
ing instabilities. One finds two phases with enhanced 
pairing correlations: 

• Dl phase: U < 2V - f3J 2 K /J H and V < 0. 

Here the dominating instability is towards conventional 
BCS singlet pairing, with correlations 

{Oss{x)Oss{x'))~\x-x'\- 1/K < . (101) 

• D2 phase: 2V-/3 J%/ J H < U< -2V, V < 0. 

In this region triplet pairing shows a power-law decay at 
large distances 

(0* TS (x)0* TS (x')) ~ \x - xT 1/K ° , (102) 

and is the dominating instability in the ground state. 
It follows that the composite singlet pairing operator 

Oc-"ss ' i d enne d i n Eq. J^SJl, also builds up large- 
distance correlations: 

ko^rg ao \x)o^t\x')) 

~ (-iy(0* TS (x)0* TS (x'))(n*(x)n*(x')) 

- (-l) e x \x~x'f 1/Kc . (103) 

We have here used Eqs. (JH3J and 1)102(1 . together 
with the property that the Neel order parameter, with 
(n(x)n(x')) = (n z (x)n z (x')) — constant, defines the 
"out-of-plane" direction z along which the triplet pair- 
ing correlations are enhanced. Similarly, for a diagonal 
stripe one would have, using Eqs. H94fl and (|102|l : 

(*)0H(s')> ~ I* " A~ 1,K ° ■ (104) 
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However, as shown by Coleman et aW°», k = odd- 
frequency pairing is likely to be intrinsically unstable, 
and hence is not expected to compete with the conven- 
tional triplet pairing mode. 

2. Away from half-filling: n e 7^ 1 



We now turn to the physically more relevant case of a 
stripe with n e ^ 1, assuming that n e is within the range 
where the stripe forms an antiphase domain wall between 
the Neel configurations (n e < 0.6 in the strong-coupling 
limit Jk ~ Jh, according to Refs. 56 and 57). 

Away from half-filling the charge sector is always in 
the (weak coupling) Luttinger liquid metallic regime, con- 
trolled by the bare value of the Luttinger liquid parame- 
ter K c . The phase diagram now splits into four sectors: 
A, B, Dl and D2 (see FIG. 5), separated by the line 
U — 2V + [3Jk/Jh = (where the spin gap closes) and 
by the crossover line U + 6V = (K c — 1) separat- 
ing the metallic phases with dominating density-density 
correlations at K c < 1 from those with dominating su- 
perconducting correlations at K c > 1. 

When U - 2V + j Jh > the spin sector flows 
to strong coupling with the excitations condensing at 
(</„) = yJn/i>K s . This kills off the CDW as well as 
the "in-plane" SDW^-correlations. The SDW 2 and TS 2 
correlations, on the other hand, survive the spin field or- 
dering, and one finds 

(O z SDW (x)O z SDW (x')) ~ \x-x'\- Kc , (105) 
(0* TS (x)0* TS (x')) ~ \x-x'\- 1/K \ (106) 

Therefore, in the 

• A phase: U - 2V + (3 J 2 K / J H > 0; K c < 1, 

the SDW 2 correlation is the dominating instability in the 
system (with TS 2 subleading), while in the 

• D2 phase : U - 2V + <3J 2 K / J H > 0; K c > 1, 

the triplet pairing (TS 2 ) fluctuations dominate (with 
SDW 2 being subleading). 

When U - 2V + (3 J\j Jh < the spin sector flows to 
the other strong coupling regime, with the spin excita- 
tions condensing at (ip' s ) = 0, the charge sector remaining 
in the weak-coupling metallic region. The correlations 
that exhibit algebraic decay are those of CDW and SS 2 
excitations, and one finds 

(Ocdw(x)Ocdw(x')) ~ \x-x'\- Kc , (107) 
(0* ss {x)0% s {x')) ~ \x-xT 1/Kc . (108) 

It follows that in the 

• B phase: U - 2V + /3J 2 K /J H < 0; K c < 1, 



the CDW correlation is the dominating instability in the 
system (with SS subleading), while in the 

• Dl phase : U - 2V + (3 J 2 K / J H < 0; K c > 1, 

the conventional singlet-pairing BCS fluctuations are the 
most dominant (with the CDW fluctuations being sub- 
leading) . 

In the A and Dl phases the composite CDW 2 and SS 
order parameters also build up large- distance correla- 
tions. In exact analogy with the D2 phase at half-filling, 
one obtains for a collinear stripe: 

(^W0f-c«(^)) - (-i)" x \ x ~ A~ Kc , 

(109) 

and 

(ot7s /ao) (z)C7 s /O0 V)> - (-i) £ x \x - X '\- 1/Kc ■ 

(110) 

For a diagonal stripe a composite CDW also develops, 
with correlations 

(C^l^f-c^^')) ~ \x - A~ K ' ■ (in) 

The fate of the zero- momentum composite singlet is more 
uncertain, however, considering its intrinsic instability 70 . 




FIG. 5: The ground state phase diagram of the stripe elec- 
tron system at n e 7^ 1. The solid line corresponds to a 
narrow metallic phase with gapless spin excitation spectrum, 
and separates two different spin gapped phases: the SDW Z 
and/or TS Z phase and the CDW and/or SS phase, respec- 
tively. The dashed line marks the crossover from a metal- 
lic phase with dominating density-density correlations into a 
phase with dominating pairing correlations. All instabilities 
shown in the phase diagram exhibit a power law decay of 
correlations. Subleading instabilities with correlations which 
decay faster then the dominating ones are indicated in brack- 
ets. As discussed in the text, the enhanced SDW Z (TS Z ) 
correlations coexist with composite CDW (SS) instabilities. 

We summarize our findings in FIG. 5. We have to 
stress again the weak-coupling nature of the phase di- 
agram. Higher-order corrections will modify the shape 
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of the phase boundaries. However, more important are 
strong coupling effects. In the case of strong values of 
the Kondo lattice interaction, one may expect additional 
phase transitions due to the finite band width of the ef- 
fective stripe model. Such effects cannot be traced within 
the continuum limit (infinite band) approach used in this 
paper, and will require numerical studies. 



D. The A phase away from half- filling: A scenario 
for nonconventional superconductivity? 

Of the phases considered, the A phase away from 
half-filling is of particular interest as the conditions 
U — 2V + PJk/Jb > 0;K C < l;n e ^ 1 are expected 
to apply to a generic stripe in a cuprate material: Most 
experiments 72 indicate that the stripes in the cuprates 
are intrinsically metallic, with no commensurability gap 
even when the stripe order is static and strong, as in 
the Nd-doped materials^. Although a precise specifica- 
tion of the coupling constants is beyond present-day tech- 
nology, the Coulomb interaction among the stripe elec- 
trons is expected to dominate other couplings (electron- 
phonon, dopant potentials, interlayer fields,...), implying 
the bound U - IV + PJ 2 K /J H > 0, with K c < 1. Un- 
fortunately, the electron dynamics on time scales where 
stripe fluctuations can be neglected (for which our model 
may apply) is still to be searched out experimentally, and 
there are as yet no " hard data" against which we can con- 
front our results. 

The A phase is dominated by a conventional SDW 2 - 
instability together with a composite CDW, coexisting 
with two subleading superconducting instabilities: con- 
ventional triplet pairing and composite singlet pairing 
(breaking parity and time reversal). This is different 
from the well-known scenario of a "spin gap proxim- 
ity effect"— where pair-hopping between a stripe and a 
spin-gapped insulating environment "infects" the stripe 
with the gap, resulting in a conventional CDW insta- 
bility, with a subleading singlet pairing channel. In the 
case where stripe fluctuations are sufficiently slow that 
they can be treated as "quasi-static", the CDW insta- 
bility can be shown to be suppressed by destructive in- 
terference between neighboring meandering stripes, leav- 
ing the singlet superconducting instability as the lead- 
ing ona^fi. The singlet order parameter on each stripe is 
then assumed to become correlated across the sample via 
inter-stripe "Josephson" coupling, leading to supercon- 
ducting long-range order below a critical temperature. In 
contrast, in our scenario singlet superconductivity (with 
the added property of breaking parity and time reversal 
symmetry) would require the suppression of the leading 
SDW 2 and (composite) CDW instabilities, in addition to 
that of triplet pairing. As we shall see below, quasi-static 
fluctuations do not perform this trick. Rather, meander- 
ing stripes living on a collinear backbone tend to phase 
lock so that (conventional) triplet pairing comes out as 
the leading effective instability. In the case of a diagonal 



structure, the (composite) singlet pairing correlations (if 
at all present; cf. our discussion after Eq. (|MJl . survive 
the slow stripe fluctuations, and coexist with the triplet 
pairing channel. Whether a complete theory — treating 
stripe fluctuations and the one-dimensional electron dy- 
namics on equal footing — would change our picture in 
favor of singlet pairing remains an open question. 

Leaving for future work the problem if and how long- 
range superconducting order may emerge from an A- 
phase type instability when stripe fluctuations are fully 
included in the analysis, there are still several issues that 
need to be addressed: 

• Is the spin gap sufficiently large for the instabilities 
in the A phase to survive at finite temperatures? 

• What happens when taking into account the fact 
that the antiferromagnetic environment as seen by 
a stripe is not that of two semi-infinite domains 
but is rather made up of two finite-width domains, 
separating the stripe from its neighbors? 

• How do transverse stripe fluctuations influence the 
A-phase instabilities? 

• What about possible long-range interactions among 
the stripe electrons? 

Let us discuss these questions in turn. 

The size of the spin gap. The A phase corre- 
sponds to a strong-coupling regime g s = — \g±\ which 
is reached after a crossover from weak coupling (where 
9s = \g±\)', cf- Sec. III. A and FIG. 3. Because of the 
crossover, the spin gap opens slowly, and it is a priori 
not obvious that it will suffice to sustain the A phase 
in presence of thermal fluctuations. To find out, we 
use <|63[) to integrate the RG equations (|62|l . and iden- 
tify the length scale where becomes of order unity. 
This scale - where the perturbation is of the same or- 
der as the fixed point Hamiltonian and renders the the- 
ory noncritical - defines the correlation length £ s of the 
electronic spin sector. Using |<7j_(^ ;j )| ~ 0(1) 3> \gos\ 
in the integrated scaling equation for g±, we obtain 
£ s = aoexp[(7r/2 — arctan(gos/<5g))/27r<5.a], where Sg = 
\J — g^ s . The associated spin gap A s = frv^ 1 is 
thus given by 



h — exp 
ao 



tt/2 - arctan(g 0s /v / .9o_L ~ 9o s 



27T V9o± ~ 9os 

(112) 

with go s = ao(U — 2V)/irv s and go± = oq{U — 2V + 
0Jjc/ Jh)/^v s , as given in Eq. JSBJ). There is consider- 
able latitude in specifying the parameters entering (|112Jl , 
but choosing v s ~ 10 5 m/s, a ~ 5A , U- 2V ~ 10 _1 eV, 
and J^/Jh ~ 10~ 4 cV — all within reasonable bounds 
— we obtain from l|112f> a spin gap A s corresponding 
to a temperature of about 500K. We conclude that the 
gap is robust, and is expected to sustain the A phase in 
the relevant temperature range (where stripe formation 
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is possible^). It is interesting to note that an estimate 
of 500K is about "right" if one were to identify the spin 
gap with the pseudogap observed in the underdoped 
metallic phase of the cuprates— . A considerable amount 
of local pairing would then be present well beyond the 
superconducting transition temperature (in this region 
determined by the onset of global phase coherence). 

Finite-width antiferromagnetic domains. A 

stripe in a real material is not embedded in a 2D an- 
tiferromagnet, but is separated from its neighbors by 
finite-width domains, or spin ladders, with a finite an- 
tiferromagnetic correlation length £,afm- In order for 
the A phase to survive in this more harsh environment 
the spin gap must develop on a length scale shorter than 
£afm (where the stripe electrons can still sample local 
magnetic order). This implies the condition 

& < Ufm , (H3) 

with £ s the spin correlation length on the stripe. The 
zero-temperature correlation length Ufm for S = 1/2 
spin ladders with an even number of legs ng has been 
calculated analytically 74 : 

Ufm « 0.5a n e a68 "' (1 - 0.73n e ), (114) 

where, as before, ao is the lattice spacing. Using that 
£ s = hv s /A s , our estimate from above, A s /ks ~ 500K, 
together with l|113(l and (|114fl . imply that the A phase 
survives for even-leg ladders with ng > 4. As suggested 
by the Monte-Carlo data inS, a four-leg ladder with Jh 
of the order of 10~ 2 eV may support the A phase up to a 
temperature of T ~ 50K. By increasing U — 2V and/or 
Jk, the spin gap grows, allowing for the A phase to per- 
sist at the lower bound ng = 4 for even higher tempera- 
tures. 

The case of an odd-leg ladder requires further analysis. 
Now, with S = 1/2, the Berry phase in l |18| l contributes 
a nonvanishing topological term to the effective action 
for the spin ladder, implying a diverging spin-correlation 
length £ — > oo but with no antiferromagnetic order over 
large distances. However, in the weak coupling regime 
the topological term is effectively inactiveZi, and as a 
consequence there is no distinction between gapless and 
gapful ladders on length scales shorter than Ufm- It 
follows that the condition 11131) is the same for even- 
leg and odd- leg ladders. Although we cannot rigorously 
exclude that nonperturbative effects may carry over to 
the stripe electrons on length scales larger than Ufm, 
it seems improbable considering the fact that the spin 
sector of the stripe develops a mass at a length scale 
which is shorter than and independent of Ufm ■ As such 
the mass is already well-established at the scale where 
non-perturbative effects from the ladder may come into 
play. 

Before concluding this discussion we wish to add 
two more notes: First, the analysis sketched here is 
strictly valid only for a collinear stripe structure. As we 



commented upon in Sec II. B, the unusual boundary con- 
dition implied by a diagonal stripe orientation influences 
the Berry phase in a nontrivial way. This may produce 
a non-negligible feed-back on the localized spins when 
considering a finite- width or spin-ladder environment, 
possibly changing the A phase in an unexpected way. 
Secondly, one should note that the composite order 
parameters for a collinear structure, defined in Eqs. I|9U|) 
and H93|) , decay faster with a spin-ladder environment as 
compared to the case of two semi-infinite 2D domains. It 
follows that at distances shorter than Ufm - where our 
construction is still expected to be valid - the algebraic 
decay of the Neel-order correlations produce a faster 
decay of the composite correlations compared to the 
conventional ones. This is different from the case where 
a one-dimensional electron gas is coupled by a Kondo 
lattice coupling to two non-interacting antiferromagnetic 
Heisenberg spin-1/2 chains (i.e. with ng = 1). As 
shown in Ref. 30, the composite order parameters (c-SS 
and c-CDW) here induce the dominant instabilities. 
This reflects the fact that the model with an ng — 1 
environment renormalizes to a fixed point different 
from ours, instead belonging to the universality class of 
chirally stabilized liquids^, with no opening of a spin 
gap. It is here important to note that if one were to 
increase our model parameters U — 2V and/or Jr- so that 
(|113fl and 1)1141) were simultaneously satisfied for ng = 1, 
our construction would break down. In particular, the 
assumption of a weak coupling RG scheme - as exploited 
in Sec. III. A - would be violated. In addition, the 
pronounced role of the Berry phase for m = 1 is likely 
to invalidate the construction of the effective spin sector 
model in l|31l) . as paths in (|17|) away from the extremum 
now enter the stage to influence the ground state also 
at short and intermediates length scales. The ng = 1 
environment is therefore not expected to be covered by 
our approach, and must instead be studied by other 
methods, such as that advocated in Ref. 30. 

The quasi-static limit: the effect of slow stripe 
fluctuations. To study the effect of slow stripe fluc- 
tuations on the order parameter correlations of phase A 
we consider a 2D array of static stripes and take an equi- 
librium average over their meanderings. Adopting the 
notation in Ref. 10, we use a coordinate system in which 
the stripe array runs along the ^-direction (collinear or 
diagonal on the lattice), with a transverse displacement 
of a stripe in the y-direction labeled by Uj{x), with j 
indexing the stripe. Introducing the arc length Lj(x), 
measuring the distance along the j th stripe to position 
x, 

Lj(x)=J dx'y/l + (d x . yj (x'))*, 

we infer from Eq. I|78|) the expression for the SDW Z order 
parameter on a meandering stripe: 

O z SDW (j,x) = O z SDW (j,x)+h.c. (115) 
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with 



O 



SDW 



ti, x) 



(116) 



The coupling of the SDW 2 to that of a neighboring stripe is of the formSi 
H S dw ~ J dxg[Ayj(x)} (o ] SDW * (j, x)6 SD w- (j + 1, x) + h.c. 



Jdxg[Ayj(x)] sin y/2n K a <p' a (J 



^2nK s if' s (j + l,x) cos ^2TrK c Aip' c (j,x) - 2k F AL J (x) (117) 
I 



where g(Ayj(x)) is an ^-dependent coupling constant. 
Neglecting possible overhangs of stripes we have here de- 
fined Ay(J, x) = y(j+\, x)—y(j, x) > (and similarly for 
A(f' c and ALj). By integrating out the stripe fluctuations 
yj (x) in powers of g one obtains an effective Hamiltonian 
of an equivalent rigid system, with a renormalized cou- 
pling 

g ~ (g[Ay 3 }) exp(-2k 2 F ((AL 3 ) 2 )) + 0(g 2 ), (118) 

where (...) denotes an average over meandering stripes. 
As argued in Ref. 10, since the signs of ALj(x) are 
randomly distributed along the distance x, one expects 
ALj(x) to grow as a random walk: ((ALj(x)) 2 ) ~ It 
follows from l|118(l that the interstripe coupling between 
SDW 2 's can be ignored in the thermodynamic limit. 

The conclusion that transverse stripe fluctuations 
cause destructive interference between SDW 2 's on neigh- 
boring stripes clearly applies to any k ^ order parame- 
ter: The exponential suppression in i|118|) can be avoided 
only if the momentum transfer multiplying ((ALj(x)) 2 } 
is identical to zero. Thus, the dephasing effect operates 
also for the A-phase composite CDW's, Eqs. 1(901) and 
(191(1 (with k — 2k F + 7r/ao [k = 2k F ] for a collinear [di- 
agonal] stripe backbone), as well as for the staggered k = 
it /a® composite singlet pairing, Eq. ((93(1 . for a collinear 
structure^!. In contrast, the conventional (k = 0) triplet 
pairing channel (|88|l survives the stripe fluctuations. 

We caution the reader that the argument, adopted 
from Ref. 10, is valid only in the quasi-static limit where 
stripe fluctuations are sufficiently slow to be treated via 
an annealed average over (static) meandering stripes. 
Moreover, the argument is strictly valid only in the ther- 
modynamic limit. The complete problem, where the 
quantum dynamics of mesoscopic stripes is treated on 
equal footing with the intrinsic Luttinger liquid instabil- 
ities remains unsolved. 

For complementary views of the physics of meandering 
and fluctuating stripes we refer the reader to Refs. 78-82. 



What about long-range electron-electron inter- 
actions? This question becomes critical when realiz- 
ing that the Coulomb interaction on an isolated stripe is 



poorly screened, given its insulating environment. Neigh- 
boring stripes may provide metallic screening over a fi- 
nite range, but our " assumption" that this range is of the 
order of a lattice spacing - implicitly built into the ex- 
tended U—V Hubbard model in © - may not be realistic. 
Still, having included a nearest-neighbor repulsion V in 
the model, we do obtain some information about the ef- 
fect of the poor screening: As can be gathered from Eqs. 

, l|105|) and (|106|) , the presence of V > enhances the 
SDW 2 instability, whereas the TP 2 gets weaker. In the 
case of an arbitrary finite screening length k s , Schulz 83 
found that the large-distance correlations are governed 
by a modified charge parameter K c ~ 1/ Vhi« s . In the A 
phase, this again gets translated into stronger (weaker) 
SDW Z (TP Z ) correlations as compared to the case with 
only local on-site interactions ~ U (cf. again Eqs. I|105|) 
and 1106(1 s ! 84 . For small J\l Jh (as assumed here) one ex- 
pects this result to apply also in the spin-gapped A phase. 
We conclude that as long as there is a finite screening 
length k s present, our results - using the U — V Hub- 
bard model - should remain at least qualitatively valid 
for length scales > k s : The A phase supports an SDW 2 
instability with a subleading TS 2 (which, however, gets 
weaker as k s grows larger). Taking into account the finite 
lengths of the stripes would introduce yet another scale 
into the problem (cf. the discussion above), requiring a 
more sophisticated analysis. 

In this context we wish to remind the reader that we 
obtained the induced spin-interaction in 1)3 1|) by a saddle- 
point approximation, with n x d T n in (|27|) locked to the 
small but fast oscillations of the £-ficld. Fluctuations 
away from the extremum are expected to produce an ef- 
fective retarded interaction in the spin sector, similar to 
what happens in the charge sector of a Luttinger liq- 
uid when integrating out electron-phonon interactions 85 . 
Whereas the local fluctuations of the Neel field could pos- 
sibly be treated in analogy with the phonon problem 85 , 
an analysis of the dispersive antiferromagnetic excita- 
tions (which may produce a long-range tail of the re- 
tarded interaction) requires a novel theoretical approach. 
We have to leave this problem as a challenge for the fu- 
ture. 
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IV. TWISTED ANTIFERROMAGNETIC 
DOMAINS 

So far we have have been concerned with the ideal 
situation of perfect phase-antiphase Neel configurations 
surrounding the stripe, giving rise to the effective U(l) 
symmetric model in (13 1 1) - In this section we generalize 
the discussion to the case where fluctuations twist the 
Neel configurations with respect to each other, breaking 
the SU(2) spin-rotational symmetry completely. 

How does the Neel order parameter change across the 
antiphase boundary defined by the stripe? Let us take 
the boundary along y = 0. The simplest situation is that 

(n(x,y)) = f(y)z (119) 

with f(y) = — /(— y). This is the situation we have con- 
sidered so far; that the rotational symmetry about the 
spin z axis is preserved, and that the Neel order parame- 
ter simply decreases across as we approach the antiphase 
boundary and is reversed on the opposite side. 

A second possibility is that the Neel order parame- 
ter rotates along (for instance) the spin x axis as it ap- 
proaches the antiphase boundary. In addition, we can 
let it rotate around the spin z axis along the boundary, 
resulting in the following form for (n(x,y)}: 

(n(x,y)) — (cos(a) sin(27rqx), cos(a) cos(27rga;), sin(a)), 

(120) 

where a = a(y) is an odd function with a(±oo) = ±7r/2. 
The parameter q is the rotational pitch along the an- 
tiphase boundary The special case q = corresponds to 
a "collinear" spin texture, and if q ^ 0, we find a phase 
that is topological, in the sense that the Neel order pa- 
rameter covers the spin space with topological density of 
Aitq per unit length of the phase boundary. 

These phases are all further subdivided by their sym- 
metry under reflections through the antiphase boundary. 
If the lattice points are arranged so that the line y = 
contains lattice points, the actual configuration is a con- 
figuration under reflections through y — and the an- 
tiphase boundary is site centered. If y = ±l/2ao con- 
tains the lattice points, the configuration is even under 
reflections through y = and the configuration is bond 
centered. 

These various stripe scenarios have all been investi- 
gated as candidates for spin textures associated with 
striped antiphase boundaries. Although the investiga- 
tions are not conclusive, it is fair to say that neither 
experimental nor theoretical investigations suggest that 
anything but q = stripe configurations should be con- 
sidered candidates for a spin texture. In fact, simple 
calculations 86 suggest that the spin texture that appears 
to be energetically favored is in fact the site centered 
collinear stripes along (10) that are odd under reflections 
through the antiphase boundary. For the diagonal (11) 
stripes, the bond centered and site centered have energies 
that are almost identical. Furthermore, the spin domains 



of these stripes are all extremely narrow, as measured by 
the width of the functions a(y). This supports our use 
of an effective ID model for a site-centered stripe. 

To study the effect of a completely broken spin- 
rotational symmetry we shall confine our attention to 
the case of a collinear stripe (along the x-direction) away 
from half-filling (cf. sec II. A), and with repulsive inter- 
actions among the electrons (U — 2V > 0). As before, 
we denote the two insulating semi-infinite domains sur- 
rounding the stripe by A and B respectively. We may 
assume that the Neel directions (x, y) and (x, y) 
are parameterized as in Eq. I|120|) . with q = and with 
a a slowly varying function of y, except across the stripe 
where a changes sign. To connect to the local mean-field 
picture used in Sec. II. A we shall simply think of two 
fixed Neel directions (nS A ') and (n^ B ' > ) which differ by 
an arbitrary angle 7, < 7 < it, across the stripe (with 
7 = 7r for perfect phase-antiphase Neel configurations). 

The construction of the effective model proceeds along 
the same lines as in Sec II. A However, since ni ^ 
—rir when 7 ^ n, we must now pay attention to pos- 
sible contributions from the coupling Hn of the lattice 
spin density to the Neel order parameters, Eq. 12111 . 
By inspection it is easy to verify that the staggered spin 
density (— l) r A r entering Eq. I|21|) vanishes in the contin- 
uum limit away from half-filling, implying that for this 
case Hn does not come into play, even when 7 =/= tt. 
Taking (n^) — > (n(a;)( A )) ~ z, as in sec II. A, it 
follows from l|27|l that the effective electron-electron in- 
teraction on the stripe mediated by the Kondo exchange 
with the A-domain remains the same as before, as given 
in Eq. (gSJ: 

H^} = J dx(:JtJl: + :J Z R J Z R : +JZJ R \+J y L J R ). 

(121) 

Turning to the i?-domain, we rotate the coordinate sys- 
tem so that its Neel direction {n^ B \x)) lies in the yz- 
plane: 

(n( B >) = sin7 y + 00S7 1, < 7 < n , (122) 

with 7 the angle w.r.t. the z-axis defining the Neel order 
direction of the A-domain. The piece of the electron spin 
density J that survives the projection onto the plane in 
which the uniform components live (cf. the discus- 
sion in sec II. A) is given by 

Jx = J (J- (n (B) ))(n (B) ) 

= sin 2 7 J y J y +cos 2 7 J Z J Z (123) 
+ cosj sinj (J y J z + J z J y ). 

Decomposing J — > J l + Jr, and using the chiral identi- 
ties 

TV jz , jz jy _ n 
J L/R J L/R ~r J L/R J L/R — U > 

Jl/r^l/r = Jl/rJl/ri ( 124 ) 
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it follows from l|27(l and l|124|) that the effective electron- change with the spins in the B domain is given by 
electron stripe interaction mediated by the Kondo ex- 



H (B) h) = - a ° J K 

n int \ !) 2J H 



J dx{:J[Jl: + :J R J R : +JfJ% + cos 2 7 J*J* + sin 2 7 J|^ -cos 7 sin 7 (J* J* + J£j£)}. (125) 



The induced interaction terms l|28|l and (|125|l are now to be added to the spin Hamiltonian H s in (|12|) (which describes 
the spin sector of the stripe electrons decoupled from the environment). Writing H s + H^t + Hint '(t) = ^(7), we 



obtain: 



ff.( 7 ) = 2™, I { :./,./,: • : J R J R : -flo,( 7 ) J£J^ - Sox-M - (<7o± - £^ sin2 7) J L J « 



a / 



K 



4irv s J H 



cos 7 sin 7 (J£jjj + J£J*)}, 



with 



W7) = — (C/-2y + sin 2 7 ^), 
and with go± and w s given in Eq. (|32|) . Introducing the rotated currents 

Jfl/Z = J fl./L, j V R,L = COS ( ) J R/L + M9)J R /L, f R /L = 'M^Jr/L + ™s(e)J R/L , 

we can write H s ("f) on diagonal form w.r.t. spin components by choosing 9 = —7/2: 



where 







= 2nv s J etej :f L f L : 




a 
ttvf 




gosh) 


a 
ttvf 


(U - 2U + sin 2 ( 7 /2) 


50/(7) 


_ a 

TTVF 


sin 2 ( 7 /2)^, 


v s 


= v F - 





(126) 



(127) 



(128) 



(129) 



T 2 

2J H ' 



Introducing auxiliary fermion fields i? M (x) and L u (x), 
connected to the rotated currents j l L / R {x) by 



?l = \'- l \o 1 uv L v :, fR = \--R\o Jl llv R v -., i = x,y,z. 



(130) 



(131) 

an application of the bosonization dictionary in Sec. 
II.C (with L M -> L n , R n -> J* /L -> j^ /L ) produces, 
as expected, a sine-Gordon model for the corresponding 
bosonic spin field <p s and its dual i? s . Using the same 
The four first terms in l|129|) are of the same form as procedure to bosonize also the last term in l|129|) — 
the spin Hamiltonian (|31|l for the ideal phase-antiphase which is new — we finally obtain 
( 7 = 7r) problem, but with a 7-dependent coupling gos( 7 )- 



(7) = T / dx {( 9 -^) 2 + ( a -<) 2 + ^-«cos(V8^7#) + -^|- K cos(v / 8^7 1 ^)}, (132) 
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where — to leading order in the coupling constants — we have: 

2(K S -1) = .g 0s ( 7 ) = ^([/- 2 y+-^sin 2 ( 7 /2)), 

2nm ( ^ = 30± ( 7 )^^([/-2y+-^(l-isin 2 ( 7 /2))), 
2™f' = 30/(7) = - — #sin 2 ( 7 /2), 

v s = v F -^(U-2V+^). (133) 



Thus, the spin sector of the stripe electrons in the pres- 
ence of twisted antiferromagnetic domains is described 
by a generalized sine- Gordon model containing a cosine 
of the dual spin field §' g . The presence of this term al- 
lows for spin-nonconserving processes, reflecting the com- 
plete breaking of the spin rotational symmetry (down to 
an Ising like Z2 x Z2 x Z2 discrete symmetry, with one 
Z2 factor for each spin component in the rotated basis 
(|128|l ). This model was first studied by Giamarchi and 
SchulaS, who used it to account for spin-orbit and elec- 
tronic dipole-dipole interactions in quasi-one-dimensional 
conductors. Its perturbative RG equations (valid in the 
limit of small parameters gos(j), 9o±(l) an d 30/(7) are 
most easily derived by exploiting the operator product 
expansion^ for the rotated currents in Eq. ^ 128(1 . using 
(|129fl . and one obtains: 

dg s 
d£ 
dg±_ 

dl 

di 



2 2 
9f~9±, 



~9s9± , 
9s9f- 



(134) 



The renormalized couplings gi = gi(j,£),(i — s,_L,/) 
are connected to the bare parameters 301(7) in l|134ll by 
9i{lA = 0) = 301(7), with I = ln(a/ao), where a is a 
renormalized length. Focusing on the case with repulsive 
electron-electron interaction (U—2V +J^/2Jh > 0, with 
K c < 1), we read off from (jl34jl that go s (-f) > 0. For this 
case the RG equations <|134[) support two strong-coupling 
massive phases: 



(i) 9os(l) < 9ox(l) ~ l3o/(7) I : 
g s -> -00, \g\_\ ->• 00, 3/^0, (135) 

(") 3o s (7) > 9o±(j) ~ l3o/(7)l : 

9s -> 00, \g f \ -> 00, g ± -> 0. (136) 



and 



By inspection, using (|134fl . we can label the two regions 
by the range of the twist angle 7: 



(i) 0< 7 <7r/2, (ii) 7r/2< 7 <7r. 



(137) 



In region (i) \g s \ and \g± \ increase upon renormalization, 
while <?/ scales to zero. The dual spin field is hence irrel- 
evant and the picture emerging is qualitatively the same 
as in the A phase of Sec. III.C.2: the SDW 2 correla- 
tion is the leading instability (critical exponent K c ), with 
TS Z being subleading (critical exponent 1/K C ). Note, 
however, that the ordering tendency is now along the 
z-direction of the rotated frame, defined by the trans- 
formation in l|128fl . It follows that for a twist angle 
< 7 < 7r/2, the leading (subleading) instability is that 
of a spin density wave (triplet pairing) oriented along 
a line in the plane spanned by (n^) and (n/ B )), and 
making an angle 7 /2 with (n^). 

Turning to region (ii), the renormalized coupling con- 
stants g s and \gf\ are seen to flow to infinity, while g± 
scales to zero. For our case, where 30/(7) > 0, it follows 
that d s gets ordered with expectation value (?S S ) = 0, 
while the spin field tp s stays disordered. By running down 
the list of bosonized order parameters in Sec III.B.l, 
we pinpoint the leading (subleading) instability as that 
of the SDW y (TS y ) correlation, with critical exponent 
K c (1/K C ). Note that also this result refers to the rotated 
frame defined by (|128fl . Translating back to the original 
frame we infer that for a twist angle ir/2 < 7 < 7r, the 
leading (subleading) instability is again that of a spin 
density wave (triplet pairing) oriented along an axis con- 
tained in the (n^) — (n^ s ')-plane, but now making an 
angle 7 /2 + tt/2 with {n^). 

Two comments are here appropriate. First note that 
the axis which defines the enhanced SDW and TS cor- 
relations makes a 7r/2-jump as the twist angle 7 passes 
through 7r/2. This may appear reminiscent of a spin-flop 
transition induced by a magnetic field, as seen in certain 
quasi-lD materials 87 . However, there is no physical ef- 
fect associated with the jump other than a sudden change 
of the direction along which the SDW and TS correla- 
tions are enhanced. Secondly, when 7 = 7r we do recover 
the result for the A phase in Sec III.C.2 (ideal phase- 
antiphase configuration). Note that the presence of the 
dual cosine-term in the effective model (|132|) simply re- 
flects the fact that we are now working in the rotated 
frame. By undoing the rotation in (|128|l we immediately 
recover the standard sine-Gordon model, as given in l|52|) 
for a perfect phase-antiphase configuration. 
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To summarize; we have shown that the opening of a 
Neel twist angle 7 it across the stripe does not remove 
the instabilities of the A phase found in Sec. III.C.2. The 
conventional leading (subleading) instability is still that 
of a spin density wave (triplet pairing), but now tilted 
w.r.t. to the two Neel- directions. 

By inspection of the composite order parameters de- 
fined in (|89|l and (|92|) it is clear that their amplitudes 
will drop as 7 decreases from its largest value 7 = n, 
defining a perfect phase-antiphase Neel configuration. It 
is here interesting to note that when 7 ^ n there is room 
for additional composite order parameters. Besides the 
trivial variations of (SSJ} and (E2J> where S (A) - S {B) 

, we can construct two composite vector order 

parameters: 

O c sdw ~ Osdw x (S (A) x S (B) ) (138) 
- SDW x (n^> x n< B >) , 

and 

O c _ts ~ Ots x (S {A) x S {B) ) (139) 
Ots x (n< A > x n^), 

where in the second lines we have dropped the small con- 
tributions from the ^-fields (which can be neglected in the 
correlation functions, cf. Sec. III.B). These composites 
have the same symmetries as the corresponding conven- 
tional order parameters, SDW and TS, respectively. In 
particular, (D c -ts is even under time and spin rever- 
sal, while odd under a parity transformation. It follows 
that for a sufficiently rigid spin texture with 7 away from 
and 7r, triplet pairing comes in two guises: one con- 
ventional mode and one composite mode, both carrying 
zero momentum. These pairing modes have enhanced 
correlations along well-defined directions, orthogonal to 
each other. As 7 — > tt (perfect antiphase boundary) 
or 7 — > ("in-phase" boundary), the composite triplet 
channel gets suppressed and one is left with the pairing 
instabilities characterizing the A phase in Sec. III. Else- 
where we shall explore the consequences of this intriguing 
possibility. 

V. SUMMARY 

In this article we have analyzed the problem of a one- 
dimensional electron liquid ( charge stripe ) embedded in 
a two-dimensional antiferromagnetic insulator, and cou- 
pled to it via a weak spin exchange. Using a mean-field 
type construction, the spin exchange gets encoded as an 
effective anisotropic spin interaction among the electrons. 
This interaction is shown to be marginally relevant un- 
der certain conditions (in a renormalization group sense) , 
producing enhanced pairing and density fluctuations in 
the electron liquid. 

For realistic values of the model parameters - assum- 
ing a screened Coulomb interaction for the electrons, and 



a conduction band away from half-filling - the dominant 
instabilities are towards a conventional spin density and 
a composite charge density wave, coexisting with sub- 
leading conventional triplet and composite singlet pair- 
ing correlations. Taking into account the slow transverse 
fluctuations of a stripe, the triplet pairing instability is 
expected to turn into the dominant one. While the mag- 
nitudes of the conventional instabilities do not change 
with the relative orientation of the Neel directions in the 
two domains surrounding the stripe, the composite cor- 
relations have largest amplitudes when the stripe forms 
an antiphase domain wall in the antiferromagnet. With 
the possible exception of the composite singlet pairing 
mode, the instabilities are found to be insensitive to the 
spatial orientation of the stripe on the underlying lattice 
(collinear or diagonal). 

Our study has been motivated by a wish to under- 
stand the role of spin exchange between stripes and their 
environment in the cuprate superconductors. As is well- 
known, these materials are exceedingly complex systems. 
It is unlikely that a simplified model like ours - where sev- 
eral aspects of the problem have been simply ignored - 
can produce accurate predictions about experimental re- 
sults. We have discussed a posteriori some of the effects 
not included in our model: the role of slow stripe fluctu- 
ations and interstripe couplings, the finite width of the 
antiferromagnetic domains encountered in real materials, 
as well as the expected poor screening of the Coulomb 
interaction along a stripe. Other aspects of the problem 
are yet to be addressed; most importantly the possible 
appearance of a retarded nonlocal spin interaction com- 
ing from large-amplitude fluctuations of the Neel order 
parameters. Also, the finite-size and boundary effects 
implied by the mesoscopic scale at which the stripes live 
need to be carefully studied. 

The virtue of our "stripped-down" model is that it al- 
lows us to carry out a well- controlled analytical study, 
and as such it could serve as a stepping stone for more 
detailed investigations. The model predicts that a stripe- 
environment spin exchange under certain conditions may 
produce instabilities towards pairing of electrons. This 
result is obtained without the assumption of a pre-existing 
spin gap in the insulating environment. When including 
effects from the fluctuations of the stripes the symme- 
try of the dominant pairing instability appears to come 
out "wrong" — as judged against available experimen- 
tal evidence from the cuprates. Does this imply that 
a stripe-environment spin exchange plays no role for su- 
perconductivity in these materials? Or, could the picture 
change when properly building more dynamical elements 
into the model? To answer these questions will require 
more work. 
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